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under Identical Quenchants, Postetch Ambients and Encapsulation 
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p- & n-Type silicon surface barrier detectors (p-SSBD & n-SSBD) have been fabricated by varying the three significant 
parameters viz., quenchants, postetch ambients and encapsulations. Postetch exposures of two days’ duration to deionized 
water for n-type and that to dust-free air for p-type have been found essential for SSBD technology. Quenchant 
experimentations lead to their ranking in decreasing order of superiority as HNO,, H,O, HF inn-SSBD and HNO,. HF. H,O 
in p-SSBD. Stain film left as a result of quenching appears to have excess of chemisorbed fluoride group. Encapsulation of 
edges and free rear surface both in p- and n-SSBD reduces the leakage current. It has been applied prior to ambient exposure in 
p-SSBD and after this exposure in n-SSBD. Leakage currents of less than | x 10 °A(n-SSBD) and 2 x 10 ~*A(p-SSBD) have 
been obtained for reverse biases up to 100 V. The sample detectors yield energy resolutions of 30 keV fwhm (”-SSBD) and 
50 keV fwhm(p-SSBD) for 5.47 MeV z-particles. Results indicate dominant dependence of SSBD technology on chemisorbed 


water and oxygen coverages at different stages of fabrication. 


1 Introduction 
Interfacial layers on silicon produced as a result of 
chemisorption are the decisive factor in SSBD 
technology. Quenchant, postetch ambients and 
encapsulation are the three fabrication parameters 
which control the composition and thickness of these 
.. layers and hence the detector response. Considerable 
investigations have been attempted on n-SSBD with 
reports suggesting a number of contradictory 
a, schedules. Successful barrier formations on p-type 
silicon have also been reported’ ~* but detailed 
scientific explanation is not attempted in them. It is felt 
that comparative evaluations of these fabrication 
parameters on p- and n-type silicon will afford a 
reliable technique to probe into the unresolved 
problems, besides leading to optimization of 
fabrication schedules of these detectors. 


L.1 Quenchant 
_ Etching and quenching leaves the first aver layer 


called stain film, the composition and thickness of 
which seem to be controlled by the quenchant. 
Deionized water has been commonly used as 
sen chant for n-SSBD but its universality is 
by Goradia* who claimed that HF and 


Dy avencings oo, produce good eto 


her's resu * indicate that this stain film consists of 
truct patios or silicon hydride and their 
ness lies in th ik gp 


ts. Still thicker films are 


superiority of quenchants in decreasing order as 
advanced by him is HF, H,0, HNO, for n-SSBD. 
Clansenn’, however, predicted this stain film to be 
consisting of some chemisorbed fluoride group. Walter 
and Boshart® correlated stain film presence to 
deterioration in n-SSBD and suggested its removal by 
surface treatments. Obviously, this reasoning would ~ 
lead to complete reversal of quenchant rankings 
advanced by Goradia®. So this controversy needs 
investigation as the results will suggest the best suited 
quenchant, besides establishing the chemistry of the 
stain film. In’ case of p-SSBD it is HNO, which is 
universally accepted as quenchant’ *, but the 
explanation for this has not been attempted in these 
reports. The relative superiorities of other quenchants 
and the influence of this stain film are yet to be 
explored for p-SSBD. 


1.2 Postetch Ambients 

These are characterized by the gases and vapours 
present in atmosphere or to which the etched 
specimens are exposed intentionally. In n-SSBD some 
reports suggest a definite improvement in the quality 
when the specimens are exposed to either water or dust 
free ait. But the controversy exists regarding the stage 
at which these exposures are to be given. Goradia‘ 
claims that only the postetch exposure to moisture is 
sufficient whereas Siffert and Coche® emphasize the 
necessity of post-evaporation exposures to air only. In 
view of this controversy and attempted 
encapsulations at different stages, the study of these 
exposures becomes necessary. In case of p-SSBD, the 
postetch exposures to atmosphere or deionized water 
are reported to be minimized'**. However, thick 
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Fig. 1—Metal-silicon surface barrier detector 


oxygen coverages on p-type silicon at this stage have 
been obtained by some workers, viz. Mathew er al.' by 
HNO, at quenching stage, Maksimov er al.’ by 
storage in discharge tube and Protsenko er al.* by 
treating it with HNO;. These coverages can also be 
obtained by exposing the surface to dry and dust-free 
air. Our attempt has been to evaluate the utility of such 
an exposure for different quenchants and different 
encapsulations. 


1.3 Encapsulation 

A major source of leakage current in SSBD is the 
surface leakage current. It flows due to the shunt path 
provided by the two space charge regions, i.e. one 
under the front metallic electrode (M,,) characterized 
by diffusion potential V, and the other under free 
surface, i.e. edges and rear side (R,.) characterized by 
surface potential V, (Fig. 1). Reports have appeared 
regarding control of this current either by 
encapsulation'® or by certain surface treatments®, It is 
claimed° that in n-SSBD if the M,.-R,. junction is 
reverse biased by keeping |V/y| larger than |V,|, then the 
difference |V,— V,| controls this leakage current. In the 
present work attempts have been made to incorporate 
both the encapsulation technique and the postetch 
exposures to suitable ambients for minimizing this 
current in m-SSBD. Reports so far do not mention 
encapsulation to p-SSBD; so investigations have been 
made to test its utility in them also. Ambient sensitivity 
of completed detectors is well known. The role of 
encapsulation in its control has also been explored. 


2 Experimental Details 

From 10 kQ-cm resistivity silicon of «111» 
orientation, round slices of surface area 3 cm? and 
thickness | mm were obtained. The same ingot was 
used to ensure identical bulk properties. Lapped slices 
were nickel! plated to obtain back contact by soldering 
tinned copper wire. The slice was cleaned in ultrasonic 


6 a 


thorough rinse in deionized water. Hence three 
different procedures were adopted to test the utility of 
encapsulation. The first (E,) consisted in encapsulating 
the entire back portion, the edges and the small 
periphery of the front of the specimen in picine. These 
portions were thus protected from etching. In the 
second technique (E;), these portions of specimen were 
exposed to etching. But encapsulation was applied to. 
p-type wafer just after etching whereas for n-type it was 
employed after postetch ambient exposure. In the third 
procedure (E;), the completed detector Le. after the 
front metallic evaporation, was mounted on perspex . 
disc by using picine at the small periphery of the rear 
side and thus most of the R,, region was kept exposed 
to laboratory ambients. The results indicated the 
superiority of E, over other procedures and so it was 
followed in the present work. 

The specimen was then etched in CP-4A etchant 
with acids of electronic grade. Quenchings were carried 
out separately for different samples by gradual 
addition of each of the quenchants HF, HNO, and 
deionized water. The specimen at this stage was 
flushed inside the etching solution by a stream of 
deionized water and was thereafter thoroughly rinsed 
in it. For studying the influence of ambients, 
identically quenched n-type specimen were subjected 
to minimum possible exposure to any ambient, one 
day exposure to dust-free air and exposures to 
deionized water for different durations. The deionized 
water was changed six hourly. The postetch exposure 
schedule selected for p-type specimen was: 
(a) minimum possible, (b) one day to deionized water 
and(c) to dry and dust-free air for different durations. 
Gold was then vacuum evaporated for front contact in 
both p- and n-type specimens. /-V characteristics were 
recorded using a Phillips GM 6020 microvoltmeter 
and employing short-shielded cables. Some of the 
detectors were tested in o-spectroscopy experimen- 
tation. Resolution curves using a thin and : 
uniforr **'Am source were plotted at 1(~? Torr a 
pressure for N, and P, detectors at 30 V reverse bias. 


3 Results a 

3.1 n-Type Silicon Surface Barrier Detectors ety “* 
Fig. 2 incorporates the reverse characteristics of ; 

these detectors quenched by different agents. Leakage 

currents below 10 -*A up to 100 V reverse bias were s 


obtained both for HNO,- and H a 
detectors. HF quenching led = poor die 
these recorded threefold larger — ae b 
besides breaking down at low bias of 50 V. } 
appears to be the best suited q ach aa 
— iority of the quenchants in decreasi 
order as HNO, H,O, HF. Fig. 3 displays the st 
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Fig. 7—Reverse characteristics of p-SSBD under different 
encapsulation procedures 


of n-type silicon. Improvements were observed with 
increasing exposures to dry air. Deterioration was 
marked when this period exceeded two days. Hence the 
optimized two days’ exposure to dry and dust-free air 
appears a suitable substitute to the other procedures of 
oxidation. Results of encapsulation experiment on p- 
SSBD are consistent with those obtained on n-SSBD; 
(Fig. 7). The first procedure (E,) did not produce a 
barrier at all. Of the other two, the second (E,) appears 
to be the best suited as it recorded leakage current less 
than2 x 10 *A upto 100 V, whereas in E; it exceeded 
1 x 10~’A even at 50 V reverse bias. Second 
procedure (E,) also enhanced the ambient insensitivity 
of the completed detector. Thus encapsulation of rear 
and edges of the sample improves performance of p- 
SSBD if it is applied just after etching. In n-SSBD it 
was employed after the postetch exposure, 

results indicate energy resolutions 
of 30 keV fwhm (N, detector) and 50 keV fwhm (P, 
detector) for 5.47 MeV %-particles at the source-to- 
detector spacing of 3 cm and at 30 V reverse bias. Thus 
performance of both p- and n-SSBD fabricated by 
optimized schedules was reasonably good. 


4 Discussion 


a Silicon Surface Barrier Detectors 
ae aa of quenchants obtained, i.e. HNOs, 
H,0, HF, is in complete disagreement with superiority 
advanced by Goradia°. It indicates that with 
increasing thickness of the stain film, leakage currents 
increase and hence deterjoration in detector quality. 
This stain film appears to have a fluoride group which, 
being electropositive induces donor levels and reduces 
barrier height on n-type silicon. Superiority of HNO; 
over other quenchants is because of its capacity to 
replace this fluoride group by initial oxygen coverages 
resulting in induction of acceptor states. The necessity 
of postetch exposure to deionized water is consistent 
with the results of Goradia®. This leads to water 
chemisorption on entire silicon surface which is known 
to induce donor levels and reduce the surface potential 
on n-type specimen. During metallic evaporation most 
of the chemisorbed water is dislodged from M,, region 
because of vacuum pull and impinging metallic: 
atoms’! resulting in an increase in V,. It cannot occur 
in R,, as this region is encapsulated and so the reduced 
V. is retained. Thus the condition |V4| larger than | V,| is 
obtained to limit the surface leakage current. Excessive 
exposure, however, may result in penetration of water 
molecules in the oxide layer of silicon thus reducing V, 
permanently. Such penetrations have been re- 
ported'*'? on germanium surfaces. The utility of 
encapsulation and superiority of E, procedure over 
others is also easily understood on the above model. 
Both in E, and E; the rear and the edges are exposed to 
deionized water but in E, these are encapsulated in 
picine prior to the exposure. So reduction of V, does 
not occur in E,. Superiority of E, over E, is because in 
the former the chemisorbed water from free surface 
cannot be dislodged during evaporation due to 
encapsulation, but in the latter the desorption of water : 
from this region does take place because of vacuum a 
pull. So the reduction of V, attained due to ambient 
schedule is retained in E, while in E; it undergoes an 
increase thus enhancing the leakage current. Fi 
the ambient insensitivity of completed detector, too, is _ 
improved in E, over E, because the reduced V, in the 
former cannot be altered by any ambient expc 
of R,. region in E, is, however, prone to changes 


displacement of water by the Oxygen atoms as most of 
R,. is not encapsulated. 
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role of thicker oxygen coverages on silicon is studied. 
Experimental data on germanium surface potential 
and the observed deterioration of n-SSBD by thick 
oxygen coverages’ indicate that the initial coverages 
enhance the overall acceptor character thereby 
reducing surface potential, and the additional 
coverages introduce donor-like states and hence 
imerease the surface potential. This reversal of 
character may be attributed to the following: (a) The 
additional oxygen coverages produce complex oxides 
with binding energies different from covalent bonding. 
It causes redistribution of surface charge and may 
result into donor levels. (b) Silicon surfaces have both 
acceptor and donor-like impurities, though the former 
is richer in p-type. The observed predominant 
diffusion of acceptor-like impurity atoms from surface 
to oxide-gas interface may leave the surface with 
relatively larger donor-like impurity. (c) Impurities 
incorporated during growth of oxides also introduce 
intermediate donor-like states'*. Hence the superiority 
of quenchants in p-SSBD is determined by its capacity 
to produce oxygen coverages. But the detector 
formatic eee ease 708s! uty 


ts that stain fiim 


ie group fees ee being quite thick in this 
© Coveney agree : 3 


increase during postetch exposure. The difference LA 

— V,| does not increase and hence the leakage current 
cannot be limited. In E, low V, due to initial oxygen 
coverages is retained because R,. is encapsulated but 
Vs in M,, increases due to additional oxygen 
coverages. |V,—V,|, being the largest in E,,yields the 
least leakage current. Insensitivity of completed 
detector is also improved in E, over E, because low V, 
is retained permanently in E, whereas in E, it changes 
due to displacement of adsorbed atoms of R,, by more | 
active ambient atoms. 
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The effect of addition of borax, nickel phosphate (NiPO,), manganese phosphate [Mn,(PO.)2], sodium dihydrogen 
orthophosphate (NaH ,PO,) and potassium arsenate on the growth characteristics and growth habits of KDP crystals at room 
temperature (25°C) has been studied. Single crystals of KDP with various concentrations of these additives have been grown by 
the Holden’s rotary crystallizer (slow evaporation) technique. The crystals obtained are of excellent optical quality and of fairly 
large size. The relative importance of the addition of growth-aiding agents like borax to KDP has been brought out. The 
presence of the additives and dopants in the grown crystals has been detected by the chemical and the X-ray analyses. 


1 Introduction 

The various physical properties of the KDP-type 
crystals and their isomorphs which fit into the 
XH ,PO, group of electrooptic materials were explored 
thoroughly by a number of workers'*. However, the 
effect of dopants and additives on the growth 
characteristics, electrooptic and photoelastic proper- 
ties of KDP crystals has not been studied in detail. 
Keeping in view that materials like borax. NiPOx,, 
Mn,(PO,),, NaH ,PO, and potassium arsenate would 
have some effect on the growth characteristics, 
electrooptic and photoelastic properties of KDP 
crystals, the present study was taken up, and the results 
obtained are presented in this paper. 


2 Growth of KDP with Borax 

Brice* suggested that borax could be used as a 
growth aid for KDP- and ADP-type crystals but did 
not give any details of the effect of borax on the growth 
rate and growth habit of KDP crystals. In order to 
Study the effect of borax on the growth rate and growth 
habit of KDP and its isomorphs, a systematic study 
was initially made on the solubilities of both ADP and 
KDP in triple distilled water with various 
concentrations of borax at different temperatures’. 

The characteristics of KDP were studied by taking 
three saturated solutions of KDP. each in 100 ml of 
triple distilled water containing 5, 10, 15g borax and 
Prepared at 36C. After cooling them to room 
temperature (25'C) they were transferred to a set of 250 
ml beakers. The temperature of the room (being a 
closed one) was fairly constant. KDP seed ee of 
known dimensions were hung in all the three solutions 
and growth was allowed to take place by slow 
evaporation technique. The increase in length of th 
crystals along [001] (along Z-axis) was monit ‘ 
periodically using a cathetometer arrangement Th 
results obtained are shown in Fig. |. The following 
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conclusions were drawn from the growth rate along 
[001] of the crystal and by visual observation. 

The growth rates for pure KDP and KDP 
containing 5 and 15 g borax in solution were almost 
same while the sample with 10 g borax had a higher 
growth rate. A small tapering was noticed in samples 
with 5 g borax whereas samples with 15 g borax did not 
show any tapering at all. In the case of samples with 
10g borax, the prism faces were observed to grow 
uniformly. The growth rates along other crystall- 
ographic axes (i.e. perpendicular to [001 ]) were high in 
the case of samples with 15 and 10g borax compared 
to those of pure KDP samples. No significant 
difference was observed in the interfacial angles with 
5 g borax and pure samples. 


GROWTH ALONG Z AXIS 
PURE KOP 
KOP WITH Sg BORAX 
KDP WITH 10g BORAX 
KDP WITH 159 BORAX 


LENGTH, cm 
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Fig. |—Growth rates of pure KDP and KDP with borax crystals 
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K DP single crystals with 5, 10 and 15 g borax with 
excellent optical quality of dimensions 45 x 12 x 12 
mm were obtained using Holden’s rotary crystallizer 
(slow evaporation) technique ina period of about 40 
days. 


3 Growth of KDP with Nickel Phosphate, Manganese 
Phosphate, Sodium Dihydrogen Orthophosphate and 
Potassium Arsenate 

Byteva® studied the effect of pH on the growth habit 
of ADP crystals in the presence of Fe** and Cr°* and 
concluded that the solubility of FeCl, in saturated 
ADP solution at a given temperature is very much 
dependent on pH. Davey and Mullin’ studied the 
growth of (101) faces of ADP crystals in the presence of 
ionic species and indicated that Fe** ions get 
incorporated in the lattice. They suggested that a 
complex of the type Fe(H,O), (OH); is formed and 
occupies interstitial sites between two rows of H,PO; 
ions parallel to the [010] direction of the (502), (701) 
and (702) faces. Torgesen and Jackson® studied the 
action of Cr? * on ADP and indicated the possibility of 
incorporation of the impurity along the growth fronts 
of the layers rather than by absorption on certain 
crystal planes. Fontanberta’, based on_ their 
Mossbauer studies, confirmed that Fe** ions do get 
into the lattice probably through the (0.25, 0.22, 0.12) 
and (0.75, 0.22, 0.12) sites. 

Davey and Mullin’ reported the growth of ADP 
crystals in the presence of manganese ions; they 
discussed the effect of manganese ions on the growth 
rates of (010) faces. They explained the kinetic effect in 
terms of absorption of the Mn** complex and the 
interaction of this with the spiral growth process. 


2.50 


LENGTH, cm 


GROWTH ALONG Z AXIS 
1 PURE KOP 

2 0.1% NICKEL PHOSPHATE 
DOPED KOP 


3 0.2% NICKEL PHOSPHAT 
DOPED KDP 


) 200 
th cates of KDP and nickel phosphate doped KDP_ 
; piney eis 


3.00 : 
GROWTH ALONG X.AXIS 
oe 1 PURE KOP 
2 0.025% MANGANESE PHOSPHATE 2 
DOPED KoOP 
250F 3 0.05% MANGANESE PHOSPHATE : 
DOPED KOP 
2.25 4 0.1% MANGANESE PHOSPHATE 


DOPED KOP 


175 


LENGTH,cm 


0.50 


0 40 80 280 


120 160 200 240 
TIME, hr 


Fig. 3 Growth rates of KDP and manganese phosphate doped 
KDP crystals 


To study the effect of nickel phosphate (NiPO,) on 
the growth characteristics of K DP, saturated solutions 
of KDP with 0.1 and 0.2% (by weight) NiPO, were 
prepared at 35 C. These solutions were allowed to cool 
down to 25 C. Seed crystals of KDP of known 
dimensions were introduced into the solutions. Then 
the crystals were allowed to grow by slow evaporation 
technique under controlled conditions. The growth 
rates of the crystals along [001 ] were monitored every 
25 hr using the cathetometer arrangement. The results 
are presented in Fig. 2. 

The growth rate along [001] was found to increase 
with the increase in concentration of NiPO, in KDP 
solution. The shapes of the pyramidal faces in the case 
of NiPO, doped samples were distorted and the 
interfacial angles decreased. It was observed that the 
tapering was severe in the case of 0.2% NiPO, doped 
KDP crystals. The tapering angle was found to 
increase with the ageing of the solution. In the case of 
0.1% NiPO, doped KDP crystals, a small tapering was 
observed, but there was no substantial increase in the 
tapering angle with time. Attempts were made to study 
the effect of higher concentrations of NiPO, on growth 
characteristics of KDP crystals. However, it was found 
difficult to dissolve more than 0.2% NiPO, in KDP 
solution (100 ml) by mere addition and agitation. 

Solutions of KDP containing 0.25, 0.05 and 0.1”, (by 
weight) manganese phosphate [Mn,(PO,),] were 
prepared. A comparative study was made on the 
growth rates of the crystals, and the results are shown 
in Fig. 3. It is evident from Fig. 3 that except in the case 
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of 0.05% Mn,(PO,),, the growth rate of KDP 
decreased with the addition of the dopant. This result 
is quite contrary to that of NiPO4 doped KDP crystals. 
The KDP crystals doped with higher than 0.1 /o 
Mn,(PO,), had higher growth rate initially. But with 
the ageing of the solution, the growth completely ceased. 
The effect of 1% (by weight) sodium dihydrogen 
orthophosphate (NaH ,PO,) on the growth characteri- 
stics of KDP crystals was studied under the same 
conditions as above, and the following conclusions 
were drawn. There was no change in the growth rate of 
1% NaH,PO, added KDP crystals compared to pure 
K DP crystals along the X-axis. There was no significant 
change in the interfacial angles and growth of the 
pyramidal faces compared to those of pure KDP 
crystals. The growth rates of KDP crystals with small 
percentages of NaH,PO, along the X- and the Y- 
directions were higher than those of pure KDP 
samples. Higher percentages of NaH,PO, in KDP 
solutions resulted in polycrystalline aggregates. 
Similarly, from the studies on 0.5% (by weight) 
potassium arsenate added KDP crystals, the following 
conclusions were drawn. The overall growth rates of 
KDP crystals with potassium arsenate were rather low 
compared to those of pure KDP crystals. No tapering 
was noticed in the growing samples. There was a small 
change in the growth habit of the pyramidal faces. 
_ To study the linear electrooptic and photoelastic 
effects in these mixed crystals, one has to have 
sufficiently large-sized and good optical-quality 
crystals. The following procedure was adopted to grow 
such crystals. 


Unsaturated solution of K DP was prepared in triple | 


distilled water at 40°C. To this solution, 0.1% (by 
weight) NiPO, supplied by Sarabhai Chemicals was 
added and the solution was stirred until the whole 
amount of nickel phosphate was dissolved. The mixed 
solution was transferred to a glass tank in Holden's 
rotary crystallizer equipment. The temperature of the 
solution was kept at 35 + 0.05°C. The seed holder 
pray KDP seed crystal was placed in the tank: 
Special care was taken in preparing seed crystals. A 
moderately large and well Sioned rar KDP 
seed crystal was selected. A basal plate was cut and 
washed with a few drops of distilled water. This plate 
was heated to the same temperature as that of the 
solution (i.e. 35°C) and hung by a thread from one of 
the branches of the seed holder. The seed holder was 
rotated at about 30 rpm. The small tapering observed 
in the growing crystal was minimized by periodical 
adding KOH pellets to the solution. : 
Excellent optical quality crystals of 0.1% (by weight) 
NiPO, doped KDP were grown (dimensions, 80 x & 
x 8 mm) over a period of about 70 days. By adopting 
the same technique, single crystals of KDP with 


NiPO,, 0.2% (by weight); NaH,PO4, 1% (by weight): 
potassium arsenate, 0.5% (by weight); and Mn,(PO4)2, 
0.02% (by weight) were grown with dimensions (100 x 
4 x 4mm),(60 x 6 x 6mm),(50 x 5 x 5mm) and (50 
x 4 x 4 mm) respectively. 
4 Detection of the Presence of Dopants in the Grown 
Crystals 

A qualitative determination of the presence of 
dopants was made by subjecting portions of the grown 
crystals to chemical analysis. By this method, the 
presence of borax. nickel, manganese and arsenate was 
confirmed in the corresponding doped crystals. The 
crystals were further subjected to X-ray analysis!° to 
study variation in the lattice parameters due to the 
presence of these ions. X-ray powder diffraction 
patterns of all the samples were obtained and 
variations in the lattice spacings as compared to pure 
K DP were observed. From the above observations, it 
was concluded that the additives had entered the host 
lattice of the grown crystals. 


5 Conclusions 
The KDP crystals grown with borax were generally 
large in size without any tapering, hard (as observed 
during grinding and polishing operations), less 
hygroscopic and exhibited very good optical qualities 
with minimal strains. The growth rate of KDP along 
[001] was very high in the presence of NiPO 4. This 
would imply that in order to obtain very large-sized 
KDP crystals, one could go in for NiPO, doping, even — 
though there are some practical problems in obtaining 
strain-free and untapered crystals of KDP in the 
presence of NiPO,. | 
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Lateral self-diffusion has been studied in indium thin films. A thin-film sample with a well defined radioactive region in the middle 
was prepared by neutron activation technique. An experimental set-up was fabricated to scan the concentration of the radioactive 
indium in the sample. A non-linear fitting procedure was employed to determine the self-diffusion coefficient from the experimental 
data. The diffusion data obtained in the temperature range 93-145 °C could be fitted to the equation D = 0.22 exp(—0.60 kT). The 
thin-film data were compared with the available data on bulk indium. A grain boundary diffusion mechanism has been concluded for 
the lateral self-diffusion in indium thin films in the temperature range studied. 


1 Introduction 

With the increasing cost of materials, especially 
some of the metals due to their limited availability, it is 
vital to develop more and more thin film devices in 
microelectronics, solar cells and other types of 
transducers. Thin metallic films are also used for 
making inter-connections, protective layers, etc. In all 
the applications, the films are subjected to electric 
fields which sometimes may produce current densities 
greater than 10* A/cm? resulting in ‘electromigration’ 
of atoms. In some cases there may exist thermal 
gradients to cause thermomigration. In either case of 
mass transport, the migration depends on_ the 
fundamental diffusion processes involved. 

In recent years, reliable data on the study of self- 
diffusion and inter-diffusion in thin films in normal 
direction have been reported using alpha particle back 
scattering', sputter ion mass spectrometers’, Auger 
electron spectroscopy?’ and radioactive tracer 
techniques*. Depth profiling by rf sputtering using 
radioactive tracers gives excellent results wherever it is 
applicable for determining the diffusion coefficient 
the normal direction in thin films. It was pointed out 
by Gupta er a/.° that data available on lateral diffusion 
measurements in thin films are meagre. Beniere e7 al.’ 
determined the lateral diffusion coefficient of copper in 
aluminium thin films by a sectioning technique using 
radioactive tracers. Scanning of the radioactive tracer 
concentration is a non-destructive technique for 
studying lateral diffusion in thin films®:*. Tai et al.'° 

| employed this method to determine the lateral self- 
iffusion in tin thin films. 

_ Dickey'! studied the self-diffusion in indium single 
‘rystals using radioactive tracers and sectioning 
technique. Self-diffusion in polycrystalline indium was 
studied by Eckert and Drickamer'?. We are reporting 
this paper the results of lateral self-diffusion in 


ned with a tracer scanning | 


films studied by neutron activation 


2 Method 

The schematic diagram of the sample used for 
obtaining the lateral self-diffusion coefficient is shown 
in Fig. l(a). The shaded regions contain radioactive 
indium (''*™[n) and the continuous stripe pattern is 
ordinary indium film. For obtaining the diffusion 
coefficients we have used an analysis technique 
somewhat similar to that used by Sun and Ohring!?. 
Fig. l(b) shows the ideal radioactive indium 
concentration profile of the sample before diffusion 
annealing when viewed through a slit of infinitely small 
width. From the fundamental diffusion laws, for the 
type of the configuration shown in Fig. 1(b), the 
radioactive concentration C(x, 4 at a distance x from 


(a) 


CONCENTRATION (ACTIVITY) 


DISTANCE 


Fig. 1—Schematic diagram of the thin film sample used for 
diffusion studies [(a), shaded region contains radioactive indium; 
continuous stripe pattern is ordinary indium film; (b), ideal 
radioactive indium concentration profile of the sample before 
diffusion annealing if scanned by a slit of infinitely small width; (c), 
shape of the profile in (b) if scanned by a slit of finite width; and (d), 
shape of the radioactive indium concentration profile after diffusion 
annealing if scanned by a slit of finite width) 
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the centre of the radioactive region of width 2h after 


annealing for ¢ sec is given by 


ee Sg es h-x )] 
C(x, = > en S x) + et TDi 


= Q and D the lateral 


Ae) 


where Co is the activity at x 
diffusion coefficient. 

To obtain the radioactive concentration profiles of 
the sample, we have placed an adjustable slit in 
between the detector and the sample mounted on a 
movable platform. Fig. I(c) shows the shape of the 
concentration profile before diffusion annealing when 
scanned by a slit of finite width. Fig. 1(d) shows the 
shape of the profile after diffusion. Because of the large 
slit width (50 um) employed in recording the tracer 
profiles, the experimentally observed concentration is 
C(x, t) corresponding to a region of finite width of the 
sample instead of C(x, t) given by Eq. (1). C(x, t) can be 
derived from the true profile C(x, 1) as 


x, )= | Pee hc thax = 
( 


<—™) 
where W is the width of the slit used in recording the 
profiles. 

When C(x, 2) is substituted from Eq. (1), Eq. (2) 
becomes 


e  Co fh Pl fhtx yee 
Ax oo ("3 a pipers oF: 
(x, A) 5 te en(; ae ei . (3) 


One can obtain D values by fitting the experimental 
data to Eq. (3) by a non-linear least squares fitting 
procedure. Statistical errors in recording the nuclear 
radiations have also been taken into account in writing 
the programme for fitting the experimental data in the 
present study. Since the slit widths used are large, the 
edge effects due to diffraction can be neglected. 


3 Experimental Details 
3.1 Sample Preparation 

Samples for diffusion studies were prepared in two 
Stages. In the first stage 99.999°, pure indium was 
evaporated in a vacuum coating unit at about 
10 ~° Torr on oxidized silicon single crystal substrates 
of 250m thickness. Indium was evaporated in the 
shape of | x 3 mm dot (0.1 yum thick) through a contact 
mask prepared by photoetching techniques. A set of 
samples thus prepared were subjected to thermal 
neutron irradiation in a neutron flux of 10!3 
neutrons/cm*-sec for one week in CIRUS reactor of 
the Bhabha Atomic Research Centre, Bombay 
Relative abundance of ''*In isotope is 4.3 . 
neutron absorption cross-section is 56 b. On neutron 
irradiation, ''*In converts to radioactive tgs 
according to 
'In+on > 4" In+7 
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Radioactive ''*™In decays to 1147y with a half-life of 


50 days 
114m y+! 14Tn +7 


1141 is a B emitter with a half-life of 72 sec 


1147p 1146p + B- 


These # particles are detected for determining 
radioactive indium concentration. 

After. neutron irradiation, in the second stage, 
continuous | pm thick indium film of dimensions 7 X 
3mm was evaporated on the irradiated dot-shaped 
samples through another contact mask at about 
10-° torr vacuum. Both the contact masks were 
prepared in such a way that on evaporation through 
the second mask, the edge of the continuous ordinary 
film coincided with the edges of the radioactive indium 
dot as shown in Fig. I(a). 


3.2 Diffusion Annealing 

Fig. 2 shows the annealing set-up for diffusion 
studies. Samples placed on a quartz sample holder 
were introduced into a quartz tube as shown in Fig. 2. 
The quartz tube was evacuated and flushed with argon 
several times and finally filled with argon up to 1/3 
atmospheric pressure. A furnace was designed 
specially for diffusion studies. A hollow metallic 
cylinder was placed inside the furnace into which the 
quartz tube containing the sample was introduced. . 
This metallic block was used for avoiding temperature 
gradient along the length of the sample. Temperature 
of the furnace can be controlled to an accuracy of 
+0.2°C by a PID (proportional, integral and 
derivative) temperature controller. Throughout the 
period of diffusion annealing the temperature of the 
furnace was monitored continuously on a strip chart 
recorder. 


3.3 Tracer Concentration Profiling 
Fig. 3 shows the tracer scanning arrangement 
for obtaining the tracer concentration profiles. A 
precision scanning arrangement was fabricated to 
move the sample laterally. The least count of the 
system shown in Fig. 3 is 1 um. For recording the — 
counts coming from a specified position of the sample, = 
an adjustable slit was placed in between the sample and 
the detector. An organic scintillator was used as a 
detector for detecting f-radiations coming from the 
radioactive indium. A single channel analyzer was 
‘no 4 


used for sorting out the pulses and for counting them. _ 
The detector has very low efficiency for detecting — 
gamma rays. When B-rays coming from radioactive _ 
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Fig. 2-Experimental set-up for diffusion annealing [1 and 10, 

thermocouples; 2, furnace coil; 3, temperature controller; 4, stainless 

steel block: 5, sample: 6, substrate: 7, sample holder: 8, quartz tube: 9, 
thermal insulator; and 1! 1, digital temperature indicator] 


Fig. 3—Schematic diagram of tracer scanning arrangement(|. single 
channel analyzer; 2, sample holder; 3, sample; 4, dial gauge: 5, 
detector: 6, slit and 7, reduction gear box) 


After diffusion annealing, samples were placed on a 

movable platform which can be moved laterally by the 

scanning arrangement. A slit of 50 um width was used 

and the activity at every 50 um distance along the 

I the sample was determined. The best fit for 

sntal data to Eq. (3), using IBM 370/155 
> the tracer concentration profiles. 
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4 Results and Discussion 
Fig. 4 shows the radioactive indium concentration 
profile of a sample annealed for 52 hr at 102°C. In Fig. 
4, the continuous curve is the theoretical curve 
obtained by fitting the experimental data to Eq. (3); 
experimental points are shown by small circles. The 
fitted values of dot width and diffusion coefficient are 
| mm and 1.45 x 10 ~°cm?/sec respectively. Diffusion 
experiments were performed at various temperatures 
in the range 93-145°C. The results are given in Table 1. 
The diffusion data obtained are found to obey the 
Arrhenius law witha single slope. The fitted values of 
the pre-exponential factor and activation energies are 
(0.22 + 0.02) cm*/sec and 0.60 + 0.01 eV respectively. 
Fig. 5 shows the comparison of our data with the 
self-diffusion data obtained by Dickey!! and Eckert 
and Drickamer'’ in bulk indium. Dickey studied 
indium self-diffusion in single crystals using 


1600 
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Activity (counts/min ) 


400 


0.4 0.8 12 1.6 2.0 2.4 
Distance (mm) 
Fig. 4—Radioactive indium concentration profile of a sample 
annealed at 102°C for 52 hr; continuous line is the fitted curve and 
circles represent experimental points 


Table 1—Experimental Data on Lateral Self-diffusior .n 
Indium Thin Films 


Temp. Time of Diffusion 
i anneal coefficient 
hr 10 ~%cm?/sec 

93 74 0.899 
102 52 1.45 
111 $2 3.72 
118 50 2.96 
126 46 5.07 
133 41 7.06 
145 33 8.88 
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radioactive tracer and sectioning techniques and of the pre-exponential factor obtainey by us isa low 
obtained an activation energy of 0.81eV. This has compared to the bulk vane (2.7 cm*/sec for 
been correlated to the vacancy diffusion mechanism. single crystal and 1.02cm /sec for polycrystalline 

indium) which is also an indication of grain boundary 


Eckert and Drickamer studied self-diffusion in 
polycrystalline indium. An activation energy of diffusion mechanism. 


0.77eV was given by them for the studies in the 


temperature range 44-145°C. We have obtained an 5 Conclusions — eee 
activation energy of 0.60 eV in the temperature range We can attribute the activation energy of 0.60 eV 


93-145°C, which is smaller compared to the data obtained for lateral self-diffusion in indium films as 
available in bulk. We could not extend the due to grain boundary diffusion. This method could be 
measurements to lower temperatures due to the employed wherever it is possible to produce a 
limitation of the method employed. radioactive isotope of convenient half-life by neutron 

Grain boundary diffusion in thin films has been irradiation. This method can well be employed for 
discussed by Gupta et a/.° and they showed that the _ lateral self-diffusion as well as inter-diffusion i in thin 


ratio of the activation energies for thin films (£,) and _ films. 
bulk (£)) is £,/E, = 0.5 ~ 0.7 for grain boundary 
diffusion mechanism. The activation energy of 0.60eV Acknowledgement 
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Large-grain (2-3 mm) polycrystalline n-InP was prepared by synthesis solute diffusion (SSD) technique. The material thus 
prepared was fully characterized by X-ray, conductivity and Hall-effect measurements. Room temperature measurements 
showed carrier concentration n, = 3 x 10?3/m, resistivity p = 3.46 x 10~* ohm-m and Hall mobility = 6.91 x 10°? m2/V-sec. 
The grain-boundary scattering effect was determined taking compensation into account. Using electrochemical technique the 
band gap was measured as 1.28eV. These results were compared with the standard sample obtained from Metals Research, 
U.K. From the temperature variation of carrier concentration and conductivity the activation energies were determined. The 
conductivity activation energy (E,) and the carrier activation energy (E,) were 9.09 x 1073 and 1.39 x 1072eV respectively. The 
temperature (7) dependence of mobility (4) was found to follow the relation 42 T* with x = 2.44. The mobility at 160 K was 


found to be 0.276 m2/V-sec. 


1 Introduction 

InP is a useful material for many electronic devices 
because of its suitable band structure and low surface 
recombination velocity. For the current optical-fibre 
communication systems, diode lasers and LEDs based 
on InP/InGaAsP are considered most suitable as 
optical light sources. Also high efficiency, stable photo- 
electrochemical solar cells have recently been 
prepared’ using single-crystal InP in which electrons 
are the minority carriers. 

The growth of InP is complicated by the high vapour 
pressure of phosphorus, 28 atm at 1070°C, the melting 
point of InP. Low temperature growth of InP was 
reported by the authors? in which the phosphorus 
vapour pressure problem and the probable in- 
corporation of impurity from the quartz crucible were 
minimized. But the quality of the crystal with respect to 
grain size, mobility, etc. was poor. Gradient freeze* and 
synthesis solute diffusion (SSD)* methods have 
recently been developed in order to grow high quality 
crystals without employing the expensive high pressure 
Czochralski liquid encapsulated technique. The SSD 
technique is fairly simple and can be employed to grow 
better quality crystal. This technique was first 
developed by Kaneko et al.*’* for the growth of GaP 
which was subsequently modified for the growth of InP 
[Ref. 6,7]. Preparation and characterization of large- 
grain (2-3 mm) polycrystalline n-InP by SSD technique 
is reported in this paper. | 


In SSD method indium (In) melt was held in a 
iperature gradient. Phosphorus (P) reacted with In 


melt on the surface to form InP which was transported 
by diffusion to the cooler end and was solidified due to 
smaller solubility. The experimental set-up is 
schematically shown in Fig. 1. A quartz crucible (1.D. 
=4mm, O.D.=6mm, length=3.5cm) containing In 
(SN) was placed in a larger diameter quartz tube (I.D. 
=8mm, O.D.=10mm, length=17.6cm) and red 
phosphorus (5N) was charged at the bottom of this 
tube. The system was evacuated to 1 x 10° * Torr and 
was sealed by fusing. The tube was subsequently put 
into a two-zone furnace whose temperature profile is 
shown in Fig. 1(b). The temperature of the phosphorus 
(Tp) was 450°C, corresponding to its vapour pressure of 
about 2 atm. The growth rate is dependent on the 
temperature difference between top of the solution (InP 
surface), i.e. T,°and the temperature of the growing 
interface, i.e. 7, and is limited by the diffusion process 
of the solute. The quality of the crystal also depends 
upon the temperature gradient. Since the diffusion co- 
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Fig. |\—{a) Experimental set-up for the SSD growth of InP and 
(b) the temperature profile of the two-zone furnace 
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efficient of P in In is low, a large temperature gradient 1s 
normally required to increase the growth rate. On the 
other hand, large temperature gradient produces poor 
quality crystals. So optimization is necessary to obtain 
the best results. In the present study, the following 
temperatures were used: 


T, =710°C, T,=910°C, T,=450°C 


and temperature gradient ~ 20°C/cm. 

The above conditions were maintained for 114 hr. 
After the growth run was completed, the furnace was 
cooled slowly (3°C/hr) up to 7,=850°C and 7, 
=650°C. Then the furnace was cooled to room 
temperature at the rate of nearly 30°C/hr. One end of 
the outer quartz tube was cut open to remove the 
crystal. 


3 Results and Discussion 

The formation of InP was verified by the analysis of 
the product using X-ray diffraction studies. No other 
phases were Observed within the sensitivity (1 %) of the 
instrument. Table 1 shows the comparison between 
calculated and standard d values for InP. 

The stoichiometry was determined by electron 
microprobe analysis taking single-crystal InP 
obtained from Metals Research, U.K. as reference. The 
results are given in Table 2. 

The results show that the stoichiometry is very close 
to that of the standard sample and there is no 
phosphorus deficiency in the grown sample since the 
observed differences are within the limit of the accuracy 
of the instrument. 
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Table 1—Comparison between Calculated and Standard d 
Values for InP 


[Radiation used: Cu" K, 2=1.54 A] 


Peak 20 d(calculated d(standard) hkl 
No. 
1 26.3 3.3885 3.39 111 
2 30.5 2.9308 2.94 200 
3 43.63 2.0745 2.08 220 
4 51.63 1.7696 1.771 311 
5 54.1 1.6951 1.693 222 
6 63.4 1.467 1.467 400 
7 69.85 1.3465 1.347 311 
8 71.95 1.3123 1.312 420 
9 80.1 1.198 1.198 422 
10 86.05 1.1298 1.129 $11,333 
il 9.6 1.0373 1.038 440 


a 


Table 2—Stoichiometry of SSD-grown InP 


In content 
Standard InP (obtained from 78.74 oe 
Metals Research, U.K.) , 
SSD-grown InP 17976 nets 
ain satin eae i CER 
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The grain size was determined optically and was 
found to be of the order of 2-3 mm. The band gap was 
measured using an electrochemical technique by 
Faktor et al®. The specimen was used as anode and 
platinum as the cathode in 1M potassium hydroxide 
solution (pH=11.5) electrolyte. The anode was 
‘lluminated and the wavelength of the incident light 
was varied between 0.4 and 0.89 um using a Jarrell-Ash 
monochromator and the resultant photovoltage was 
measured. A photovoltage versus wavelength plot was 
obtained (Fig. 2) and from the cut-off wavelength the 
band gap, as determined using the relation E,=1.24/4, 
was found to be 1.28eV at room temperature. 

From the thermoelectric power measurements it was 
found that the material was n-type which was also 
verified by Hall-effect measurements. Conductivity and 
Hall mobility were determined from _ Hall-effect 
measurements using van der Pauw technique. Ohmic 
contact was made by evaporating In-Sn alloy and 
subsequent annealing at 400°C for 1 min in argon 
atmosphere. The results are summarized in Table 3, 
which also shows the comparison between prepared 
and standard InP samples at room temperature. It is 
evident from Table 3 that there is a good agreement 
between the properties of the prepared sample and 
those of the standard sample grown by LEC. The 
higher mobility of the single crystal is due to the 
absence of grain-boundary scattering. For the prepared 
sample, the grain-boundary scattering effects were 
calculated taking compensation into account as 
described by Jensen’. Table 4 summarizes the results. 

From the temperature variation (160-300 K) of 
conductivity, carrier concentration and mobility, 
activation energies and temperature dependence of 
mobility were found out. Temperature (7) dependence 
of mobility (1) followed the relation xoc T * where the 
value of x was found out from log T versus log p plot. 
The results are summarized in Table 5. 


Voltage (mv) 


i 


Fig. 2— Variation of photovoltage with wa ar ie j 
termi iets velength for 1 


Table 2 aa between the Properties of Prepared and Standard InP Bi dainiss at Room Temperature 


Prepared sample Standard sample (LEC) 
(SSD) from Metals Research, UK 
X-ray diffraction d values close to 


the standard values 
and no other phases 


Grain size 2-3 mm Single crystal 
Band gap* (eV) 1.28 1.29 
: Resistivity (ohm-m) 3.46 x 1074 2x 107* 
i. | Mobility (m?/V-sec) 6.91 x 107? 13x 107! 
: Carrier cone. (n,), m~? 3 1073 2 x 1073 


*Also measured by electrochemical technique. 
S ‘Table 4—Compensation and Grain Boundary Parameters of comparatively low mobility value (6.91 x 1072 m2/V- 


= Se SSD-grown InP sec) is mainly due to grain- boundary scattering while 


ion of charged donors, Nj (m~*) 3.5x10-?3__ the effect of compensation is small (r = 1.333). ae 
ation of charged acceptors, N, (m-*) — 0.5 10? tee 
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Specific conductivity of magnesium valerate (soap) solutions in binary liquid mixtures h 
C is independent of benzene, chlorobenzene and p-xylene concentrations in methanol 
on: log A= A + Blog C where A and Bare constants and C is the 
h increase in temperature and methanol concentration in the 
nd independent. The soap is a weak electrolyte which dissociates below the CMC. 
_dissociation constant (K), activation parameters of dissociation, AH”, AG” and 


soap aggregates into micelle and the CM 


and the temperature. The conductivity data fit well in the equati 
concentration of soap in mol dm °. Constant A increases wit 


mixture but the values of — B are fou 
Molecular conductivity at infinite dilution, ¢ , 


6 November 1981; revised received 8 February 1983 


as been studied. It is found that the 


AS° have been evaluated and the effect of composition of binary liquid mixtures has been discussed. 


1 Introduction 

Conductivity measurements are the basis for some of 
the earliest theoretical developments in the field of 
micelle formation and structure. Two basic character- 
istics which distinguish surfactants from ordinary 
solute are their surface activity and their ability to form 
micelle by an association-dissociation equilibrium. 
The relationships between these properties and their 
uses in greases’ *, paints and varnishes°, catalytic 
actions®’’, etc. have been studied extensively. 

Magnesium valerate is found to aggregate into 
micelle in water and the critical micelle concentration 
(CMC), dissociation constant, K, molecular con- 
ductivity at infinite dilution, 4,, and heat of 
dissociation, AH° have been evaluated®. The 
conductivity of these soaps in alcohols has been 
investigated”"'® which confirms that CMC _ is 
unaffected by change in temperature. 

In the present paper, we report the conductivity of 
Magnesium valerate in binary liquid mixtures of 
methanol with benzene, chlorobenzene and p-xylene at 


different temperatures in order to study its micellar 
behaviour. 


2 Experimental Details 

Methanol, benzene, chlorobenzene, p-xylene and 
magnesium carbonate (Merck or B.D.H.) were used. n- 
valeric acid, obtained from Sigma Chem. Co was 
purified by distillation under reduced pressure and the 
purity was confirmed by b.p. measurements. 

The calculated amount of Magnesium carbonate was 
weighed and suspended in water. The suspension was 
heated up to 80°C and the required amount of valeric 
acid was added in portions under constant stirrin 
After evolution of carbon dioxide had ceased, the 
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magnesium valerate was purified by recrystallization 
from alcohol and then dried in vacuum. 

The procedure of conductivity measurement was the 
same as reported earlier'’. Mixtures of known 
composition were prepared by volume, taking 
precautions as far as possible to keep the solutions out 
of contact with the atmosphere. Lower percentages of 
solvent mixtures were used which allowed the distinct 
measurement of conductivity of solutions. 


3 Results and Discussion 

Specific conductivity (k) soap solutions increases 
with increase. in concentration of soap. The plots 
(Fig. 1) of k vs C indicate that a point is reached at 
which k increases at a slower rate. The concentration at 
which the inflection occurs, i.e. 0.01 M, shows the CMC 
of soap. This is independent of benzene, chlorobenzene 
and p-xylene concentrations in methanol, temperature 
and the composition of the mixture. The CMC of 
magnesium valerate, 0.01M in methanol is also 
found’? to be independent of temperature. _ 

Molecular conductivity (A) obtained from the 
measured specific conductivity of solutions decreases 


Cx 103 
Fig.1—Plot of specific conductivity k (ohm! cm~') ys 
concentration C (in mol dm~*) of magnesium valerate in methanol 
(B) in presence of chlorobenzene (A) at 35°C 
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Fig. 2—Plot of 4 vs C'? of magnesium valerate in methanol (B) in 
presence of benzene (A) at 35°C 
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Fig. 3—Plot of log A vs log C of magnesium 
: golame.s in methanol (B) in presence of p- 
= a (A) at 35°C 


the . The plots (Fig. 2) of A vs 

a’ ie tess due to the weak 
aviour of soap. 

— of soap_ solutions at 

peratures and compositions of eS oa 
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Fig. 4—Plot of C*4? vs 4“! of magnesium valerate in methanol(B) in 
presence of chlorobenzene (A) at 50 C 


dissociation of soap can be derived as follows: 


(CsHpCOO),Mg=Mg?* +2C,HoCOO7 aseha) 
C(1 —«) Ca 2Ca 

Ca(2Ca)?  _ 4C7a3 
phic P abgeee Ns gl) 


Since & is very small, ionic concentrations would be 
low; therefore interionic effects can be neglected. On 
substituting %=A/A,, in Eq.(3) we obtain the 


expression: 
3 2 
eo —- | (4) 


The applicability of Eq.(4) has been tested 
graphically (Fig. 4) by plotting C?A? vs 4" '. The plots 
are linear below CMC and the equation fails above this. 
concentration. The values of A,, and K have been 
evaluated from the slopes and intercepts of these plots. 

~ Molecular conductivity at infinite dilution 1 increases 
as the temperature increases in different solvent = 
mixtures and decreases with increase in ‘methanol — 
concentration at different temperatures. A compa = 
of A,, values of soap in different solvent mixtt cs of ae 
varying Le Soa shows that = are 2 the 


Tad. 


‘Ge tw 
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CAG? 
Table 1—Values of AH®, AG®, AS° (at 35°C), ( aT ) and 


dAS° ; in P of 
aT of Magnesium Valerate in Methanol in Presence 
6 

Different Liquids 


IS ae) e& hee 
—AG° =A? 2 oT or 


I~! kKJK ~!mol™! kJK ' kJK 
mol ' 


Vol. % of 
liquid 
present —AH® 


kJ mol”! kJ mo 


Benzene ° 
10 47.86 
20 43.07 
30 38.29 


: Chlorobenzene 


10 67.01 
76.58 
95.72 


0.18 
0.20 
0.20 


0.28 
0.27 
0.35 


O17 = 
0.20 
0.27 
0.35 


33.50. 
43.07 
57.43 
PFE»: 


se ve 1 that the tahies of —AH” vary aanilatle with 
trations: of additives as that of K. 
energy AS has been pesaet from the 


_ 


é_ 
1 
} 


~2 mol? 
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Values of entropy of dissociation (AS: °) hay re als 
been obtained from AH® anc AG’. Both free energy: 
entropy of dissociation increase with increa 
concentration of chlorobenzene and p-xylene i 
mixture, whereas in benzene it shows the re 
behaviour. The temperature coefficients of free ener 


and entropy of dissociation have been calculated ar : ae 


are recorded in Table 1. 
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The two-state distortion approximation in the impact parameter formulation is employed to calculate the cross-sections for 


He*(1s) + H(ls) + He*(1s) + H(n=2) collision process for the incident energy ranging from 16 to 2500 keV. The present 
results are in reasonably good agreement with those of the existing four-state calculation and the experimental findings. 


1 Introduction | 
a ~ In a previous investigation’ we have studied the 
excitation of the projectile He* ions in collision with 
hydrogen atoms using an approach based on the two- 
state distortion approximation? in the impact 
parameter formulation. On comparison with the 
many-state close-coupling calculation, we have found 
that the distortion calculation’ is in reasonably good 
[Se t specially in the intermediate- and high- 
region. The low computing cost and the 
| numer simplicity makes the distortion calculation 
attractive as compared to the -coupled-state 


e mentioned energy region. In the present 


the tote tal cross-sections for the Process: 
(1s Pee ee 


lation for describing collisions of systems having 
nal iru ure for both the target and the projectile _ 


we propose to apply the same ‘method! to. 


with 
aes ae <®,,| V;|®,> eg VinmeXPLi€mnS/?'] 
Enn = Em— En 
Sa SStut 
2 
V; — eee Stee ies (3) 


YA, rB, Vi2 
®,(r4., rp, =O" (r4 )O# (rg) 
n 1? =) 1 A, n B, 


the symbols being the same as used in our previous 


paper!.. . 
For 2s and 2p4 state excitations, all the Van 


involved are even in s, and on simplification one © 


obtains the distortion probability as 


fd. eee 


* ‘Pe tie Coca: mn 10 te: 
und state to the mth stat ¢ Is given 
[ caste 


ee 


1? : 7 4 
Pe ict ae 


probability 
at EP : Lu Wet acl sd at 
a ale 1 ke Bas a aot ie 
» 4 = - 
- Fa 
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ross-sections, 77; and 2» 


ae : 
(in units of a6) in the Reaction He* (1s) + H(1s) ~ He (1s) 


Table 1—Distortion Approximation C 


+ H(n=2) 
First Born 
sti Four-state 
5 Full two-state § 
bi 98 ati . mation 
qoegtn mi Phere calculation of calculation of geen 
He" et calculation Flannery Flannery 
(keV) am o =v ay 02s F2p 02s O2p 
O25 G2» 2s 2p 

16 1.04( — 2) 3.39 —1) 1.01(—2) 9.20{ — 2) 1.22(—1) 1.04 —1) pat es 
5 7.04 —2) 9.16(—1) 5.86( — 2) 4.02(—1) 4.30( —1) 3.78( — 1) sp fe 
36 1.93( —1) 1.68(0) 1.65( — 1) 1.19(0) 1.02(0) 9.03(—1) ae ee 
64 4.28( — 1) 3.35(0) 3.90( — 1) 3.30(0) 2.08(0) eget soy pt 

= 2.12(0) 3: : : 
100 5.16(—1) 4.63(0) 4.86( —1) 4.57(0) beets 
- .35(0) 4.18(0) 7.02( —1) .09(0) 

200 4.51(-—1) 4.80(0) 4.42(-—1) 4.74(0) | 

400 Be —1) 3.67(0) 2.94( — 1) 3.65(0) 6.24 —1) 3.46(0) 3.69( — 1) oe 
900 1.55(— 1) 2.31(0) 1.54(—1) 2.24(0) 2.38(—1) 2.20(0) 1.70( — 1) ee » 
2500 6.03( — 2) 1.07(0) 6.03( — 2) 1.02(0) 7.32(—2) 1.02(0) 6.25( — 2) 1.00(0) 


Note: The numbers in the parenthesis denote the powers of ten by wh 


ich the numbers are to be multiplied. 


them with the two-state and the four-state impact 
parameter calculations of Flannery’. It is found from 
Table 1 that at high energies of the projectile the 
present calculation for 2p excited state coincides with 
the first Born, the two-state, and the four-state 
calculations. On the other hand, the cross-sections for 
the 2s excited state at high energies obtained by 
applying the four-state calculations” differ from the 
first-Born, the two-state, and the present distortion 
calculations. It may, however, be noted that the two- 
state calculation of Flannery” does not differ much 
from the present distortion calculation above 200 keV. 
Since the present distortion calcuiation neglects the 
effect of back coupling, the close agreement between 
the full two-state and the present calculation indicates 
that the effect of back coupling is not very much 
appreciable for the s-state excitations. 

In Fig. 1, our results for the excitation of hydrogen 
atom by He”* ion impact in the 2s state are compared 
with the theoretical results? obtained in the first Borh, 


) 


2 
fe) 


CROSS-SECTION O5- (in units of a 


re) ' 2 3 4 S 
VELOCITY OF INCIDENT He’ ION (au) 


Fig. 1 —Cross-sections for 1s 2s excitation of hydrogen atom by 
He* ion impact in the reaction He* (1s) + H(1s) + He*(1s) + H(2s) 


the two-state and the four-state calculations in the Seton anne Peat ~-— first Born 
i annery’); — —, A 
impact parameter treatment. The full two-state results calculation (Flannery?) oe. Foch fea aul parameter 


of Flannery agree well with the present distortion 
calculation throughout the energy range of the 
projectile. The cross-sections calculated by the four- 
state method? agree reasonably well in shape with the 
present calculation; however, the magnitude of the 
four-state calculation is somewhat larger in the 
intermediate-energy region. With the increase of 
energy, the cross-sections obtained by the first Born 
approximation approach closely to the present 
calculated results. However, the four-state results 
predict somewhat larger values of cross-sections than 
the first Born and the present calculation. This 
overestimation, according to Flannery?, is due to the 
strong 2p-2s dipole coupling which has been 
incorporated in the four-state calculation. 
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calculation (Flannery?)] 


In Fig. 2, we present our results for the excitation of 
hydrogen atom in the 2p state by the impact of He * j 
and compare the results with the first Born, the two- 
state, and the four-state calculations of Flannery?. The 
two-state results almost coincide with the present 
distortion results except at very low energy side and 
this again indicates that the effect of back coupling has 
negligible effect in the cross-sections. The cross- 
sections obtained by the present two-state distortion 
calculation agree well both in shape and magnitude 
with those from the four-state calculation. At high 
impact energies of the projectile the cross-sec 
obtained by using the first Born, the two-state ar 
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CROSS-SECTION >, (in units of ag) 


0 1 2 3 4 5 


INCIDENT VELOCITY OF He’ ION (a.u.) 
Fig. 2—Same as Fig. 1, but for ls > 2p excitation of hydrogen atom 
by Be ion ‘impact in the reaction He* (1s) + H(1s) > He *(1s) + 
mio Hp) 


four-state calculations. are identical sath those 


_ obtained i in the present calculation. 

_ In Fig. 3, we compare the total cross-section seat 
ie Process (1) with: the recent experimental 
s of Aldager al.*, Young et al.°, and McKee 

1 the available theoretical predictions” 
ental results of McKee et al.!° for the 2s 
C ; have been added with the reported 
tio measurement of Young et al.? for the 
ion of the 2p state of hydrogen atom to allow 
riso r 1 for the process (1) with other total cross- 
sults. sonoma results as well as those of 
alculation pee well with the 


+1 Datta, frei nt 5401) Ss 


ui 
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TOTAL CROSS-SECTION (10 ’cm?) 


10 20 50 100 200 
ENERGY OF INCIDENT He* ION (keV) 

Fig. 3—Total cross-sections for the reaction He*(Is) + H(ls) > 
He *(1s) + H(n=2) (Theory: ——, present distortion calculation: — 
—, first Born calculation [Flannery’); ..., full two-state impact 
parameter calculation (Flannery’):; —-—, four-state impact 
parameter calculation (Flannery); —--—, first-order Glauber 
calculation (Franco*); ©—-—, symmetrized first-order Glauber 
calculation (Franco*): ---, Vainshtein-Presnyakov-Sobelman 
approximation calculation (Theodosiou®). Experiment: jg and A, 
Aldag et al.’; O, McKee eral.'°: 2, Younger al.?; @, Aldag et al.*) 


approximation are found to be in good agreement 
with the data in the energy range 25-175 keV. The 
V: inshtein-Presnyakov-Sobelman results® are also in 
good agreement with the experimental values over the 
energy range 50-200keV. Above 50keV, all the 


°. existing theories predict cross-sections with varying 


agreement among themselves. The present distortion 


calculation is, however, in neeboniahy. ee —_ = 


tie four-state aaleulaviont: =a ge Son 
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Using nonperturbative solutions, the binding energies for different 
n with r has also been studied. 


potentials. The variation of the wavefunctio 
earlier calculations. 


1 Introduction 

It is well known that there are some processes which 
although occurring outside the nucleus are dominated 
by distances which are small on the atomic scale 
(Compton wavelength distances). In atomic photo- 
effect, for example, Pratt and Tseng’ have argued that 
for a wide range of photon energies, electron Compton- 
wavelength distances are of primary importance. In 
such a case, the knowledge of screened wavefunction is 
desirable for including the effects of screening, at least 
for small r. Further, internal conversion”, threshold 
pair production® and single quantuin annihilation* 
phenomena are also characterized by small distances 
on the atomic scale. The normalization screening 
theory has recently been applied successfully to explain 
the anomalously large photo-defect cross-section in 
molecular hydrogen’. There are many atomic and 
nuclear processes which are characterized by the 
behaviour of an electron wavefunction at the origin. 
One such phenomenon is the orbital electron capture 
because, in this process only the region of overlap 
between electron and nuclear wavefunction is 
involved. 
The screened Coulomb potential also finds 
importance in atomic phenomena involving electronic 
transitions. It has been treated numerically and 
analytically by various workers using different 
methods, such as WKBJ method®, the quantum defect 
method and different types of perturbation methods. 
Realizing the utility of the screened Coulomb potential, 
several variational calculations®~!2 have been done 
and tables of energy as a function of screening 
parameters have been compiled. McEnnan et al}? 
used analytic perturbation theory for nonrelativistic 
case, while Greene and Aldrich'® applied nonper- 
ho approach to the problem. Later, Mehta and 
i Sorell oye 

; & energies for different 

atoms. 
_ When a classical charged two-particle system is 
influenced by a plasma sea, the Coulomb potential is 
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atoms have been evaluated for screened Coulomb 
The results obtained are in excellent agreement with 


replaced by a static screened Coulomb potential, the so 
called Debye potential’® to explain the interaction. 
Ray and Ray’® obtained s-matrix, discrete energies 
and wavefunctions for the s-states for exponential 
cosine screened Coulomb (ECSC) potential in Ecker- 
Weizel (EW) approximation. Recently, Dutt et al.'7"'* 
obtained bound s-state energies of an electron in 
ECSC potential by analytical method using EW 
approximation. They further proposed an extension of 
the EW approximation to treat the non-zero angular 
momentum bound-states of a class of screened 
Coulomb potentials and obtained the discrete energies 
for Yukawa potential. 

Motivated by the growing importance of the 
screened Coulomb potential, we have considered in 
this paper a more general potential of the form 


Ze? 


MO= AT + Sor) 


- 
where do is a screening parameter. 

The energy eigenvalues and the wavefunctions for 
different values of n for potential (1) have been obtained 
by an approximate and nonperturbative approach. 

On comparing our results with the corresponding 
experimental values and the results of Mehta and 
Patil'*, it is found that very good fits for bound state 
energies are obtained for n= 1 in potential (1). It is also 
seen that for this value of n, the bound state energies for 
lower values of Z have poor agreement with the 
experimental values but the agreement improves for 
the medium and higher values of Z. Also the shapes of 


wavefunctions for n=1 tally well with those of 
McEnnan et al.'? 


2 Calculation of 
The radial Schrodinger equation with potential (1) 
can be written as 
Pee See See 
| IP dr dr Al+or)* oF 


. 


r)= ER(r) 


ae 


; In Eq. (2), we have used atomic units (unit of length 
_—- %q =h*/me? and units of energy = —me*/h? and of 
=(%o)"do]) is a dimensionless quantity. 

- Setting Y(r) = rR(r), Eq. (2) is transformed to 


| -345 oe = Ey(r) ... (3) 
In Eq. (3), 


>a may also be written in the form 
~Zér~! 
= 1+0r’)In( + or") 


ry Sha O71 + 1)r2"~ 2 6 ) 
+ Fini + or")[(1 ae Tse ar 


as Vir) tends to Vir) for or” <1, i.e. 


(6r")? <1, and hence it is neglected. an vlan 

Thus Eq. (3) with the help of Eq. (5) takes es form 
a car 

| ~2a- G+6r)ind tor) 


| 
| 


q = 571+ 1)r2"~ 2 an : 
Din + SPN + OP"? = (1+ 57 j en =BHe9 


ces .. (6) 
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ges or ~ 
a M4 , ge a (7) 
‘ is “he : 
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In Eq. (9) terms involving 5? (for the region 6? < 1) 
have been neglected. Further to observe the usefulness 
of the results with n # 1 in potential (1) we restrict the 
study for those processes which are determined in the 


region r ~ |. Hence one gets the following for n = 2 
from Eq. (8) 


V1 — y)p"(y) — 03 — 2y)e'(y) 

we SE Z 

+59 ()) — 2aye'(y) + 757) 

—S(1 +9 +)7)9(y)=0 ts (10) 


and the following for n = —1 


yl — y)e"(y) + y?@'(y) — 3 + 2a)e(y) 


* 2¢(y) a (F —2a Jou 


—2ay(1 + yy) =0 PRS. . (11) 


2.1 Solutions for the Case, 7=1_ : 
Assuming ¢(\’) to be of the form 


a= Cyl = ices w+ (12) 


and substituting it in Eq. (9), one finds on equating 


coefficients of like powers of y that o = | for s-states and oe 


the C;’s satisfying the following recursion relation 


bs te 2a 2Z/6\C, ee 3) 


Cae CIGD) 


If the series (13) i is to terminate after a aalidale ’ 


pete ee ie ~ af i seodeadiees 


(MIM 4 2ay + D=226=0 oO aie ‘ _ 


with the 
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yields the following for the binding energy 
Zi 8 
eobeed 


2.2 Solutions for the Case, n =2 


Again for n = 2 in potential (1), C;'s satisfy the 
relation 
= (di+2a+ 1/2) -(Z/20) +(a/2)]2Ci San 
cnet (2i +1) 
with o=0. 


This yields for s-states and for 
M =), do =Z/6 


and 
M=1, for a, =(4 -2) 
= 1, for a, = ee 


2.3 Solutions for the Case, n= —1 
Similarly for s-states and for n = —1ando = 0, the 
recursion relation 1s 


[di+2a+2)—2(Z/d) + 2a]C; 


Le ft9) 


a (+1)(i+2) Se 
Eq. (20) yields 

ap =Z/6 

and 


(25~,) 
a, =|=;-- 
0 aa 

The above values of ay finally give the following 
expression for the wavefunctions 


vse Mt) 


M 
— »—4yor" a 
Wudlr)=e ahi, CyLIn(1 + 6r")]” * sin, Le) 
The wavefunctions given by (22) are approximate 
and their resemblance with exact wavefunctions 


so On the extent to which Ver) approximates 
VPN). 


3 Results and Conclusion 

The binding energy values (in keV) obtained from 
Eq. (17) for n=1 for different atoms with Z ranging 
between 34 and 84, and the values of the screening 
parameters chosen and the corresponding experi- 
mental values are given in Table |. 

It can be seen from Table | that the a 
calculated and observed values is poor for ines 
of Z but it improves for higher values of Z. 

In Table 2 the binding energy values calculated for n 
= —I are given. They show a reasonably 
agreement with the energy values obtained by lyer'®, 
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It may also be mentioned that the screened Coulomb 
potential (1) changes to a form considered by Mehta et 
al’* for n = —1. The energy values obtained in this 
work are comparable with those obtained by them. 

It may be noted that the method of investigation 
used in this paper involves minimum number of 
parameters and is quite straightforward sare 
compared to the methods used by other workers” *’. 
The potential considered by us has a general form and 
the comparison of shapes obtained for the 
wavefunction Wr) for n = 1 (Figs 1 and 2) reveals 
that within the interior of the atom, i.e. dr < 1, shapes 
of all the screened Coulomb potentials are similar and 
are Coulombic in nature. This enhances the chances of 
possible use of the wavefunctions Wr) as the trial 
wavefunctions. Further, there is a shift of maxima for n 
= 2 towards origin for the increasing values of Z. For 
the same n, the magnitude of the wavefunction falls off 
rapidly with increase in r for higher values of Z (Fig. 3). 
Finally we conclude with the remark that fairly 


Table 1— Bound State Energies E (in keV) with n=1 as a 
Function of Z and 6 
x é E 
Present study Experi- 
with n=1 in mental 
. potential (1) (Ref. 14) 
34 0.09 — 1.52 — 1.27 
39 0.095 — 1.84 — 1.70 
44 0.100 —2.19 —2.21 
49 0.120 —2.93 —2.79 
54 0.130 — 3.50 — 3.46 
59 0.135 —4.12 — 3.98 
64 0.16 —5$.11 — 5.02 
69 0.18 — 6.20 — 5.94 
74 0.20 — 7.39 —6.95 
84 0.22 —9.22 —9.3] 


en 
Table 2— Bound State Energies E (in keV) withn = —lasa 
Function of Z and 6 


& é E 
er a 
Present study By dispersion By pertur- 
for n= —1 in relation bative 
potential (1) (Ref. 14) approach 
(Ref. 19) 

34 0.09 — 1.52 = 1.30 —1.36 

39 0.095 — 1.84 —1.74 —1.55 

44 0.100 —2.19 -2.25 —2.12 

49 0.120 —2.92 + 2,82 —2.80 

54 0.130 —3.49 — 3.45 —3,53 

59 0.135 =—4.11 -—4.15 —4.44 

64 0.160 —5.10 —4,92 —5.22 

69 =—s «0.180 -581 — «$9 ; + ae 

74 0.200 -7.37 — 6.66 —_ 

0.220 9.20, 8.66 ’ Sey = eee 
i 
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: r 
Fig. |—Unnormalized radial function Wyir) for 1s state of 
aluminium (Z = 13) for n = 1 (Distances are in electron Compton 
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ae 


satisfactory results of binding energies for potential (1) 
are obtained only for n = | and that the utility of 
results for the cases with n + 1 is restricted for the 
processes determined in the region r =~ 1. However, 
such a restriction does not exist for the case n=1. 


References 


1 Pratt R H & Tseng H K, Phys Rev A (USA), § (1972) 1063. 


2 Band I M & Sliv L A & Trzhaskovskaya M B, Institute of 
Technical Physics, USA Report No. 262 (1970); 
unpublished order No 1314, 21/IV. 


3 Tseng H K & Pratt R H, Phys Rev A (USA), 4 (1971) 1835, 6 
(1973) 2049, 


4 Tseng H K & Pratt R H, Phys Rev A (USA), 7 (1973) 1423. 

5 Cooper J W, Phys Rev A (USA), 9 (1974), 2236. 

6 Schiff L 1, Quantum mechanics (McGraw Hill, New York), 1968. 
7 Seaton M J, Mon Not R Astron Soc (GB), 11 (1958) 8504. 

8 Lam C S & Varshni Y P, Phys Rev A (USA), 4 (1971) 1875. 
9 Roussel K M & O’connel R F, Phys Rev A (USA), 9 (1974) $52. 
10 Greene R L & Aldrich C, Phys Rev A (USA), 14 (1976) 2363. 


nH Rogers FJ, Graboske H C& Harwood D J, Phys Rev A (USA), % 
(1970) 1557. 


a] Harris G M, Phys re 125 (1962) 1131. 


238 McEnnan J, Kissel L & Pratt R H, Sad Reo bile 134979 


Indian Journal of Pure & Applied Physics 
Vol 21, July 1983, pp. 408-412 


Random Electrodynamics of Nonlinear System: 
Part I—Quartic Anharmonic Oscillator 


S SACHIDANANDAM} &I VV RAGHAVACHARYULU*t 


Bhabha Atomic Research Centre, Trombay, Bombay 400085 


Received 26 April 1982; accepted 19 April 1983 


A successful extension of the classical techniques of Random Electrodyn 


amics to nonlinear microsystems is still not 


obtained in the literature. A beginning is made in this direction in this paper. We study the quartic anharmonic oscillator as an 
illustrative example. By extending one of the approximation methods employed in the study of deterministic nonlinear systems 
to stochastic nonlinear systems, properties quite close to those given by the quantum mechanical description are obtained. The 
results partly dispel the doubts raised by Claverie and others in the validity of Random Electrodynamics in the description of 


nonlinear microsystems. 


1 Introduction 
Interest in Random Electrodynamics (RED) stems 
from the possibility that it may emerge as a classical 
covering theory for quantum mechanics. It postulates 
the existence of a force due to radiation damping 
accounting for the radiation losses and a fluctuating 
zero-point electromagnetic field. Incorporating these 
two forces in the equations of motion of any classical 
electrodynamical system, a Langevin type of equation 
of motion’ of the system in RED is first obtained. One 
can then study the solutions of these equations directly 
or develop a Fokker-Planck type of equation for the 
distribution function as the problem is stochastic in 
nature. These two approaches produce equivalent 
results for simple linear systems, essentially in 
agreement with those obtained on the basis of 
quantum theory. As regards the linear oscillator, the 
quantum results are not qualitatively different from the 
classical ones: The resonance frequency is independent 
of the driving fields and remains the same as the 
classical characteristic frequency; only the mean square 
displacements and velocities get quantized. The 
quantum oscillator is thus a classical oscillator on 
which theradiation forces of damping, and the zero- 
point fields act independently. However, the Fokker- 
Planck approach when applied to nonlinear systems 
has yielded results which are in disagreement with the 
results obtained when the system is treated quantum 
mechanically. The disagreements range from a sign 
difference in the second order correction to the energy 
of a slightly anharmonic oscillator to the catastrophic 
self-ionizing solution of the hydrogen atom which is, 
on the basis of quantum mechanics, obviously 
unphysical’. Unlike in the case of linear systems, the 
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situation is drastically different with nonlinear 

oscillators. The greater the deviation of the system 

from linearity, the more profound is the effect of the 

driving fields on the system dynamics, necessitating in 

principle, a nonperturbative approach. It is quite well 

known that such methods are not available in 

methematical literature. But if RED is to be the cover 

theory for quantum mechanics, then obviously it 
should yield essentially the quantum results. 


We hold that the discrepancy in the results is due to 
an inadequate treatment of Nonlinear Random 
Electrodynamical System from the point of view of 
mathematics rather than lack of soundness in the basic 
physics itself. 

In view of the nonavailability of exact treatments in 
the literature, the purpose of this note is to develop, at 
least in the perturbative framework, suitable 
mathematical techniques that enable us to study 
nonlinear RED systems using the Langevin type of 
equations and illustrate the method as applied to a 
quartic anharmonic oscillator. The other much more 
physical problems like that of the hydrogen atom will _ 
be considered separately. é 
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2 Preliminaries of Random Electrodynamics 

The physical idea in RED is as follows: 

An electromagnetic system, subject to the additional 
forces of RED, evolves from an arbitrary initial state, tc on 
an equilibrium state in which the electromagnetic 
energy radiated away by the system is balanced on an | 
average by the energy picked up from the ; Sse§rve ir of 
random fields. That is, under equilibrium, the r; d tio ; 


ae 


~~ 


SACHIDANANDAM & RAGHAVACHARYULU: RED OF NONLINEAR SYSTEM PARTI 


2.1 Radiation reaction force 

The damping force due to radiation reaction can be 
considered in a variety of ways: (1) in terms of self- 
action®; (2) through the canonical formalism in 
classical mechanics with the inclusion of interaction 
with radiation; (3) via the nonrelativistic quantum 
electrodynamics in the Heisenberg picture*; and (4) by 
considerations of the energy loss through radiation’. 
However, the energy loss argument is the most suitable 
one for our present purpose, since it throws light on the 
limitations of the resulting form of the force and avoids 
the idea of self-action, by leaving the source of the 
reaction force unspecified. 

Now consider a point charge in accelerated motion. 
The non-relativistic Larmor formula for the rate of 
energy loss by the accelerated charge is given by: 

dE ; 
“dr 30" 


We can picture the above loss as arising due to an 
equivalent radiation reaction force Fz. This force is to 
take energy from the particle at the same rate as given 
by Eq.(1) on an average, when we set 


i Qe? F 
— Kdt=>-3 rd 
\ F,.rdt xe |, r.rdt 


2 r Tt 
—se4|F | oh r rat ae tz) 
0 0 


We note that if the motion is periodic or nearly so, the 
first term on the right can be neglected, and we get 
2 

=i was. 3) 
In view of the neglect of the first term, it would not be 
correct to use this formula in situations involving 
arbitrary aperiodic motions of the particle. Moreover, 
even in periodic motions, the term containing r should 
be negligible compared to the other force terms so that 
the term with the third order derivative can be reduced 
to a lower order one using suitable proportionality 
‘constants’. By this procedure we suppress the 
unphysical runaway and preacceleration solutions. It 
should be emphasized that one must always bear the 
above limitations® of the Lorentz damping force in 
mind. Otherwise one would get the type of counter- 
intuitive solutions discussed in detail by Jayaratnan - 


»veff) 


F,= 


2.2 Zero-point Electromagnetic Fields (ZPF) 

As is well known, quantum electrodynamics (QED) 
leads to the concept of a vacuum state in which each of 
the modes (k./) of the free electromagnetic field has a 
zero-point energy value 1/2ha,; By virtue of the 
resulting divergent energy density. the zero-point 
energy Is eliminated in the ‘first’ renormalization. 


q However, it has been found useful, in certain contexts, 


to visualize an equivalent real classical fluctuating 
electromagnetic radiation, satisfying the homogeneous 
Maxwell equations. For example, the spontaneous 
emission by an excited atom is physically pictured as 
emission stimulated by the ZPF. In the calculation of 
radiative corrections such as Weisskopf's*® cn the self- 
mass of electron and Welton’s’ on the Lamb shift, the 
ZPF can be assumed to induce a jiggling and a 
translational Brownian kind of motion respectively, on 
the electron. Since the heuristic calculation of Welton, 
various arguments have been advanced for taking the 
ZPF as real. One should of course always keep in view 
the limitations when using these fields, the most 
important one being the imposition of a suitable upper 
cutoff to eliminate the divergences. 

A simple model for these fields with their infinite 
number of degress of freedom is that of a Gaussian 
random process with the spectral density 
2h 
S(W) == lol’ .. (4) 
3¢ 
A convenient representation, for ease of application in 
the Langevin type of differential equations, is, in terms 
of random phases, as follows: 


; = 
E(r,t)=1/2 (rs 3 e(k./) ro hd 


A=1 
exp[ —i(wt—k.r)]+c.c} aE (| 


where a,,=exp(i0,,) with 0,, being a random variable 
uniformly distributed in the interval 0 to 20, 
independently for each wavevector k, and polarization 
index 4. The usual orthogonality and closure 
properties hold good for the polarization vectors e(k./). 
They ,are: 

e(k./) .e(k.2) = 93; 


e(k.4).k=0 
and 
F vik. Ae (k.2) =) — kik? . 6) 


Further, in view of the independence of the phases for 
each mode (k./), we have the following ensemble 


average <)> over phases: 

Cay My 4) =9= (ar ds > 

and 

Kaya” Y= ok —k)0;, .. (7) 


2.3 Langevin Type of Equations in RED 

Consider an electromagnetic system which is 
subjected to the two forces of damping and radiation 
fluctuation fields besides any other usual deterministic 
fields that are already acting on it. Then we obtain a 
generalized equation in RED of the given system. The 


409 


INDIAN J PURE & APPL PHYS, VOL 


are very simil angevin 
equations so obtained are very similar to the Lang 


equations. | 

In the original Langevin approac 
motion, the damping force arises as the s : 
of the fluctuating influence of the reservoir on the 
Brownian particle. On the other hand, in the quantum 
mechanical derivation of the fluctuation-dissipation 
theorem by Callen and Welton’®, a system losing 
energy by dissipation, becomes subject to a fluctuating 
force; in particular, a charged harmonic oscillator with 
a radiation damping term given by Eq. (3), 1s acted on 
by a fluctuating electromagnetic field precisely with the 
spectrum of the zero-point field plus that given by the 
Planck formula for blackbody radiation. 

In RED, we assume classical dynamics, deduce the 
damping terms from considerations like those of Section 
2.1 and assume the existence of the zero-point fields 
independently. The resulting Langevin type of 
equation has several complexities not present in the 
original Langevin equation, not the least of which is the 
nonwhite nature of the spectrum defined in Eq. (4) of 
the fluctuating forces. 


h to Brownian 
ystematic part 


3 The Quartic Anharmonic Oscillator.in RED 

Though a general basis can be developed in treating 
stochastic nenlinear systems making use of the 
Langevin type of differential equation within the 
framework of perturbative mathematical theory, we 
illustrate the general principle by making use of the 
quartic anharmonic oscillator as an example. 


3.1 Equation of motion 
A particle moving in one dimension in a potential 
given by 


V(x)=m(w6x?/2 + Bx*/4) .. (8) 
satisfies the familiar equation 
X¥ + mex + Bx? =0 ... (9) 


Further, we assume that the particle has a charge e 
which couples it to the radiation reaction and the zero- 
point fields. Hence, incorporating the two inde- 
pendently acting fields also in the equation of motion 
defined by Eq. (9), we obtain the stochastic differential 
equation of the charged particle given by 
¥ + 6x + Bx? — t%'=e/m E° .. (10) 
where t=2/3e?/mce°, and EQ is as given by Ea. (5). 

The Eq. (10) cannot be solved exactly. However, in 
having recourse to approximation methods, one 
should keep in view the peculiar effects of non-linearity. 
The equation is very similar to that of the Duffing 
oscillator'' except that the harmonic driving forces in 
the latter are replaced by radiation fields having the 
spectrum defined by Eq. (4). 

From the literature on the Duffing oscillator, we 
know that in an iterative or perturbative calculation, 
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the ‘natural’ frequency of the oscillation has to be 
renormalized at every order of the calculation. The 
reason for this is that the frequency becomes a function 
of the amplitude and the force term. A simple way to 
take the major part of this feature non-perturbatively 
into account is through an optimal linearizational! In 
this method, the original equation is replaced by a 
linear equation with the frequency as a parameter. The 
resulting solution is used to evaluate the time average 
of the square of the difference between the two 
equations. The parameter is then fixed by requiring this 
quantity to be a minimum. We extend the optimal 
method to the Langevin type of equations by replacing 
the averages over time by those over the ensembles of 
random fields in the steady state. 


3.2 Solution of the Linear Equation Motion 
In view of the foregoing, we first solve, instead of the 


Eq.(10), the linear equation: 


A) 
where Q is to be chosen to minimize the mean square 
difference between Eq.(10) and Eq.(11). With E® as 
given by Eq. (5), the steady state solution of Eq. (11) will 
be 


¥—1¥+07?x=e/mE? 


x en x c| (12) 

0? — w’*—-itw? mn 
where we have used the dipole approximation in using 
the fields evaluated at x=0. Making use of the 
expectation values in Eq.(7) for averages over phases, 
we readily find that 


1 e2 a 


ay 


A=1 


._ hw 
fe ex(k.A)55 


: 2 
*(@ =o?) +208 


.. (13) 


Using Eq. (6) for the summation over the polarization 


index e(k,4), performing the angular integration, and 
using w? =ck?, we get 


<x?) =(h/mn)(2/3)(e2/me3) 


| * wide 
s (Q2 ie «w?)? + t2@® 
The integral is evaluated by making use of the fact that 


the integrand is sharply peaked at @=Qand replacin; 
all o's [except in terms involving (Q— @)] by Qto yield 


ve (1S) 


w. (14) 


<x?) =h/(2mQ) 
Similarly, we get 
<x)? =(h/2m)Q 


r 
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provided we eliminate the logarithmic divergence by 
assuming a suitable cutoff limit which incidentally 
provides a radiative correction to the right hand side of 
Eq. (16). 

The nth moment of x? can be similarly evaluated by 
contracting the a’(k,4), a(k’,/') factors into all 
possible pairs so that the 2n-fold integration in the k- 
space becomes a product of n simple integrations. Since 
there are (2n)!/(n!2") such pairs, we have 
(x?")> = {(2n)!/n!2"\<x2>" 


This completes the analysis of the harmonic 
oscillations of a charged particle in RED, at zero 
temperature. And the analysis can be extended to 
higher temperatures by modifying the spectral density 
given in Eq. (4) by multiplying it by coth (hw/2kT) so 
that for example <x*> becomes 
«x?>7= <x?)  coth(AQ/2kT) 


rhe, 5 Wg 


. (18) 


The above analysis of the linear oscillator in RED was 
first obtained by Marshall' and most of the subsequent 
papers and review articles'* '* in the field of RED 
have revolved around it without any essential advance 
in the understanding of the scope of RED with respect 
to realistic microsystems. 


3.3 Parameter Choice (Optimal) and Evaluation of 

Ground State Energy of Oscillator 

In order to obtain the ground state solution of the 
linear system dealt within the previous section 
maximally close to the true solution of the quartic 
anharmonic oscillator, we shall choose Q to minimize 
the quantity 


a? = (we x + Bx? —Q?x)?> 
= wa—7)?<x?> + 2B(wh -—Q?)Kx*) 


+ B?<x®> .. (19) 
On making use of Eq. (17), this becomes 
h ry 
eS: = 0>2)3 ee 
ees: - er; aaa - (52a) 
... (20) 


2f(_h \ 
“Sd Aa : 
In terms of the variable € =Q/awp , and the parameter 


y=(3/2)(hB/mw,), the function to be minimized 
assumes the form: 


zy Sy? 


2m \ 1 
se)= (fans) HO -2U-2yt+5 tat 35 rag ee 
Setting Af/0E =0 gives, 
: BE? — 26+ — &? —4yé —5y =0 2D) 


Now the roots of the equation must lie close to those 


: for the case y=0, ie. near €7=1 and ¢?= —1/3. The 


ly relevant solution must, therefore, be of the form ¢ 


eet 5. where 6 is a small quantity of order ). 
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Substituting in Fq. (22) and ingluding terms up to 
second order in »(6), we get 


3267 +(8 ~45')0 —(4) + Sy?) =0 
so that 
0 = Wo (1 + 7/2 — y7/8) 


. (23) 


(24) 
Taking the second derivative of Eq.(21) it is easily 
verified that the value obtained is indeed a minimum. 
We can now compute the mean energy of the quartic 
oscillator at the absolute zero temperature as 
CE> = 1/2m¢X?> + 1/2niwe <x?) + | 4mp <x*> 
= 1/2m<x?>(Q? +02 +3 2<x?)>) oe 4.) 
Making use of Eq. (24) and the fact that (Eo) =1/2hao 


which is the mean zero-point energy of the 
unperturbed oscillator, we have 
(E> =(Ep)(1+9/4-7°8 
3 hp 18/ hp \? 
=< E,>s 1+-—, - —|—> 2 
‘ 5 8 mg zal nes a 


The above vaiue for the mean energy in RED is to be 
compared with the value obtained in quantum 
mechanical perturbation theory given by 


3 hp 21/ ph 2 
8mw, 64\ma2 


and the Stochastic Electrodynamics (SED) value of 


3 Bh 141/ Bh \? 
— ee 
8 mus * 556 (25) Cm 


obtained by Pesquera and Claverie'?. 

An evaluation of the absorption spectrum of the 
oscillator can be performed only after elucidating the 
stationary states of the oscillator in equilibrium with 
the Planck fields in addition to the ZPF. This would 
require the solution of transcendental equations 
instead of Eq. (22). However, the renormalization effect 
of the Planck fields on the optimal frequency will be 
negligible because of the smallness of the value ol 7 (i.e. 
w2 > 3B<x?>). We can thus conclude that at 7=0, the 
first absorption peak for an applied small signal will 
occur at 
@=W (1 +7/2) w. (29) 
to the first order in 7. This agrees with the quantum 
mechanical value and has again to be contrasted with 
the value Pesquera and Claverie'* had obtained 
through the linear response theory, i.e. 


4 Discussion | 

The disagreement between Eq. (27) and Eq. (28) in 
the second order term even in sign is taken by Pesquera 
and Claverie to be one indication of the failure of SED. 
Now in quantum mechanics, as is well known, the first 
order term in Eq.(27) is the expectation of the 


(Egg (E a) 4 eee #4) 


a= (E> 1 te 


vs (30) 
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in the unperturbed eigenstate 


of the linear oscillator while the second-order one 
results from the first-order correction to the 
wavefunction. Thus, it would appear that Pesquera 
and Claverie’s method does not identify the correct 
quantum state of the system. On the other hand, that 
Eq.(26) is quite close to the quantum value is an 
indication that our generalized optimal determination 
of the frequency takes us quite close to the proper 
ground state of the system. 

Pesquera and Claverie estimate the frequency as a 
function of the energy of undriven quartic oscillator for 
use in their Fokker-Planck type of equation. This 
parallels Rauscher’s method of iteration! ° which starts 
with the solution of the autonomous system as the first 
step. The agreement of their results for the energy at 
least up to the first order shows that the system is 
basically linear. The greater the deviation from 
linearity, the larger will be the error in the approach, as 
the atomic example shows :In the case of the hydrogen 
atom there is no zeroeth approximation in which the 

- orbital frequency is independent of the orbit size. The 

____ radiation forces must hence be of decisive importance 
in determining the mean frequency of the orbital . 

_ motion. Planck’s constant enters the frequency in the 

form je*/h> whereas in the slightly anharmonic 

oscillator it appears in a small additive correction term: 


eg + = w+ 3ph/4 moe 


perturbing Hamiltonian 


| 1€ atomic problem, the mean frequency and the _ 
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when all its possible states of stationary motior 
not been elucidated by studying the equilibriun 
the Planck radiation and the zero-point fields. The 
of detailed radiation balance in nonlinear syste! 
noted by several workers’® is not likely to persist é 
such a study. We are encouraged in this attitude 
another development, though it is not related to 
nonlinear question. Several fundamental qua 
concepts such as spin | 2 and identical pa 
correlations. lying somewhat outside the mainstre 
the developments in RED reported in the literat 
far, have been shown to be related to the basic i le 
RED, in a preliminary report by Sachidanar 
Also, work on the hydrogen atom is in progres 
will be reported separately. ; : i 


References a ie 
1 Marshall T W. Proc R Soc London (GB), 284 (196 
2 Claverie P, Pesquera L & Soto F, Phys Lett A | 

(1980) 113.0 7 hag 
3 Jackson J D, Classical electrodynamics (Wi 
Edn, 1975, 786. Ws 
4 Milonni P W & Smith W A. Phys Rer A(US 
5 Heitler W. Quantum theory of radiation. 
Edn, 1954. 26. eg 
6 Landau L D & Lifshitz E M. Classical 1 
Press. London), 1962. 231-5. 
7 Jayaratnam C E, Rev Mod Phys (| 


ce 


elt 


Indian Journal of Pure & Applied Physics 
Vol 21, July 1983, pp. 413-415 


7 

. 

) Evaluation of Spectroscopic Parameters for the Pr?+ 

Ion in a Laser Liquid 

S V J LAKSHMAN* & S BUDDHUDU 

eee Laboratories, Department of Physics, § V University, Tirupati 517 502 

Received 4 May 1982 

a 

Radiative lifetimes are predicted for the fluorescent *P,,*Po and 'D, levels of Pr?* in the laser liquid SeOC1,. Quantities 
such as squared reduced matrix elements, electric and magnetic dipole linestrengths, total radiative transition probabilities and 
; radiative relaxation rates which are required for lifetime calculations are also reported. 

i 

: 


1 Introduction The right hand side matrix elements of 
Experimental and certain theoretical investigations (f*SL||U*||f%SL’) are taken from the tables reported 
of the absorption spectra of a few lanthanides in by Nielson and Koster® and the values of the 6j symbol 
SeOCl, at room temperature are available in the {JJ'/ 
Siswakure!°. Since Feats lifetimes for the )L’ ao (for 2 = 2,4 and 6) are taken from the tables 
fluorescent levels *P,,*Po and 'D, of the Pr>* ionin reported by Rotenberg e: ul.’ 
SeOCl, laser liquid are not available, the authors took 
_ up calculation of these and the procedure followed and (II) INTERMEDIATE COUPLING ; 
_ the results obtained are aches in the present is ows The calculation of the reduced matrix elements in an 
“Ss Se intermediate field is illustrated for a typical transition 
Se detest) Tntensitics Hitirsdis eae hee like 'D, > — 3F, belonging to 4/?(Pr°*) configuration. [ee 
eh cetimates of pneeaiiies for the electronic The energy matrix for J = 2 is a3 x 3 matrix sce = ae 
pee tte Pr Bein keane asm = ak SS 
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¢3p, | U4||3P) in LS coupling are easily determined 
from Eq. (2). From these and the eigenvectors (C,, Ca, 
6. Cy. C,.:and C,) obtained by solving energy 
matrices’? the reduced matrix elements 
(p;|| U4||3F5) for the intermediate coupling case are 
calculated using Eq. (4). The reduced matrix elements 
are thus transformed from the LS basis states to the 
physical coupling scheme prior to being squared and 
substituted in Eq. (1). The squared values of Eq. (4) 
were substituted in Eq. (1) and using fexpt for f.4 the 
values of 7, were evaluated by the least squares fit 
method. With these 7, (A = 2, 4 and 6) values, 
oscillator strengths f° were calculated for the 
different transitions. 


2.2 Method of Electric (S,,) and Magnetic (S,,,) Dipole Linestrengths 
Theoretical spectral intensity (f© ”) can be 
evaluated by the relation 


: iz 
itn 2 Sa 1S | (5) 


ft -370+1) On 


where n=refractive index of the medium, S,, =electric 
dipole linestrength, S,,4= magnetic dipole linestrength, 
m=mass of an electron, c=velocity of light, v 
=wavenumber of the band in cm ', h=Planck’s 
constant,e=electron charge in C,andJ =value of the 
initial state J. 


The value of S., then is given by: 
2 =) 
S.4=e| 1.08 x 10! ued 


(+1) 9 TAM WU Wa’? ... (6) 


The reduced matrix elements in Eq. (6) are 
transformed from the LS basis states to the physical 
coupling scheme prior to being squared. The value of 
Sma is then given by: 


e*h? 


h 
Sind = Ténmicz WIL + 2S ||Ws')? os AQ 


Values of ||L + 2S\|? are as given below (assuming the 
selection rule to be AJ =0 + 1): 


J’=J 
(SLI||L + 2S||SLJ’) = gJU + 1)(2J + 1) 


1 IV +1) + SS + 1)— LL +1) 
where g AOU od ae 


J’ =J-1 
(SLI\|L + 2S\\SL, J —1) 
=(S+L+J + INS+L+1-JV+S-LU+L~5) 
ET EW + LS) 
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J'=J +1 
(SLIWL + 2S||SL, J+ 1) 


_ | _s4L4+s+2(StJ+1-—) 
4J +1) 


(L+J+1—S\(S+L—J)] ot I@ 


The matrix elements calculated from Eqs (8)-(10) 
were transformed into the intermediate coupling 
scheme before computation of the magnetic dipole 
contribution represented by Eq. (7). 


3 Results and Discussion 


The absorption spectrum of Pr** in SeOClI, was 
measured by Heller’. Using his results, eigenvalues and 
eigenvectors were evaluated by solving the energy 
matrices of 4/2(Pr**) given by Dieke® and Spedding” 
The eigenvectors thus obtained were used in the 
calculation of squared reduced matrix elements || U||? 
between the electronic excited states 7P,, *Po, 'D, and 
the next lower lying levels of Pr** in SeOCI,. S,., and 
Sma Values were evaluated from Egs (6) and (7) 
respectively. Oscillator strengths were also calculated 
from the Judd-Ofelt and electric and magnetic dipole 
linestrengths methods. 


Dieke® has reported that 7P,, *Po, 'G4 and °F; 
states would exhibit fluorescence in Pr*>* in LaCl, 
single crystals. Carnall et al.'° have pointed out that 
*P,.°Po and 'D, states of Pr** exhibit fluorescence in 
LaF, solution. Since the present work is on solution 
spectra, radiative lifetimes for *P,,*Po and 'D, excited 
States have been theoretically evaluated. 


Carnall et al.'' suggest that the total radiative 
transition probabilities (A) can be calculated using the 
Judd-Ofelt intensity parameters (7;) evaluated in 
absorption measurements. The authors have followed 
the same procedure in the calculation of A for the 
excited states °P,,°Po and 'D, of Pr** in laser liquids. 


The total radiative transition probability A(wJ; wJ’) 
was calculated using the relation: 


4,,3 
Ad: WI’) = oe reer 
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The total radiative relaxation rate (A) was 
evaluated using the formula: 


AW) = ¥ AW: WJ’) 
vr 


a (12) 
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Table 1—Total Radiative Transition Probabilities (A), 

Radiative Relaxation Rates (A,) and Predicted Radiative 

veuscsmedbe! for the Levels of Pr** in SeCOCI, + SnCl, and 
SeOCl, + SbCl, 


Transition 
SLI 


7, 


A(s~') 
SeOCl,+SnCl, SeOCl, + SbCl, 
0.01 0.02 
11.81 9.42 
1104.89 301.66 
3947.33 1978.05 
23016.25 4501.97 
12543.97 8554.72 
6193.53 5208.98 
17590.60 16665.70 
8238.41 3929.60 


7.2646 4.5150 
13.77 22.15 


Value of 4,(s~')x 10>¢ 
Value of t,p(s) 


21.53 
1875.31 512.63 
4034.62 2116.98 | 

0 0 

34573.33 22261.96 

9777.77 8139.83 
0° 0 
23075.97 10951.72 


7.3288: 4.3940 
13.64 22.76 


b722 


ION 


where the sum runs over all //J’ lower in energy than 
WJ. The radiative lifetimes of a state is given as: 


tr(W) =[A(W)]-! 


Total transition probabilities A(W;: WJ’), total 
radiative relaxation rates (A) and predicted lifetimes 
(tk) for the levels of Pr*?* in SeOCI, are presented in 
Table 1. It is interesting to note that the radiative 
lifetimes of fluorescent levels are much higher (about 
two times in ?P,,>Po and five times in 'D,)in the SbCl, 
complex than in the SnCl, complex. 


. (13) 
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Schottky barrier MIS diodes have been fabricated for two 
different values of interfacial layer thickness, by a chemical method 
with Ag, Cu, Ni and Pd as the barrier metals. The barrier heights of 
these contacts have been determined and utilising these data and the 
theory developed by the authors earlier, the neutral level (~,) and the 
charge present in the interfacial layer (Q)) of such diodes have been 
calculated. 


In a recent paper’, the authors have indicated a 
method for determining (i) the amoung of fixed charge 
(Q,,) in the interfacial layer of a MIS tunnel diode and 
(ii) the neutral level (@ .) at the surface of the 
semiconductor with the help of the equations: 


Q,. 
= E ~ 4 = OF 
P n01 g Xs qDs, : (1) 
and 
Ona = E, + 1, — Go — Sox te) 
: . qD 5, 
Q,. 
qDs, 


Where Pmoi And Ymo2 are the points of intersections of 
the theoretical line a= @m— 7s with the plots of the 
experimental values of @ zg, versus (, of number of MIS 
contacts. There will be two points of intersection 
between the experimental and theoretical lines since 
there are two experimental curves, each corresponding 
to a particular value of interfacial layer thickness (0) 
and surface state density (Ds). In the above equations, 
E, and Xs represent respectively the bandgap and 
electron affinity of the semiconductor and Ds, and Ds, 
the densities of states at the surface of the 
semiconductor for two different layers of oxide 
thickness (6). Since the values of Ds, and Ds) can be 
easily calculated from the slopes of the Ms, versus 7) 
plot, it follows that Qo, and @» can be determined from 
Eqs (1) and (2), provided E, and 7s are known and ~ 
and Pmor are experimentally measured. 7 
Based on the above principle, an attempt has been 
made in this laboratory to determine the values of Q 
and p» experimentally for an epitaxial Silicon surface 
on which several metals were deposited 


416 


chemically. For the purpose of this experiment four 
metals, viz. Ag, Cu, Ni and Pd were selected for 
deposition on n-type epitaxial silicon wafers having an 
epi-layer thickness and resistivity of 20 4m and 4-6 
ohm-cm respectively. The wafer surface was etched 
and cleaned in the usual way”. Oxide layers were then 
grown on the epi-layer by two different methods, one 
by boiling the wafer in conc. HNO; for 20 min 
following the technique already reported’* and the 
other by heating the sample in an atmosphere of dry 
oxygen at a temperature of 450°C for 45 min. The 
oxide thicknesses were measured by an ellipsometer 
and found to be about 25 and 60 A respectively. The 
four metals metnioned above were then deposited 
selectively using the chemical reduction and electroless 
process described earlier?”. 
The barrier heights of the contacts thus formed were 
measured by both /-V and C-V techniques. The mean 
values of these measurements were plotted against @m 
(Fig.1), curve a being for 6 ~ 25 A and curve b for 0 = 60 
A. The work function values® were Ag(4.31 eV), 
Ni(4.74 eV), Cu(4.52 eV) and Pd(4.95 eV). On the same 
plot, the theoretical curve obeying the relation 9, =, 
— 7,(assuming 7, to be 4.05 eV) was also drawn. From 
the intersection of the latter curve with the curves a and 
b, M9, ANC P,,, were Obtained as 4.78 eV and 4.77 eV 
respectively. From the experimental plot of @,, versus 
?,, (curve a and curve b), D,, and D., were found to be 
1.26x10'?  states/em~7eV and 1.13x10!3 
0.8 


-> 
© 0.6 
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wm 
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NOTES 


- states/em ~ eV respectively. Substituting these values 
in Eqs.(1) and (2), Qo. and @ were obtained as 1.0 
; ee 10"? charges cm ~? and 0.31 eV respectively. 

__ The above method thus appears to be a simple one 

_ for measuring the charge in the oxide layer of a MIS 

tunnel diode. An examination of the experimentally 

tained curves (Fig.l) also reveals that below a 

‘icular value of the metal work-function (¢,, = 4.68 

‘the barrier height for larger values of interfacial 

layer thickness is more compared to that for smaller 

+ thickness. This trend is reversed for work-function 
i above this ¢,, value. 
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Aluminium antimonide thin films of different thicknesses were 
prepared using co-evaporation method. The variations, with film 
thickness, of the Hall coefficient Ry, Hall mobility “x, and carrier 
concentration 7, have been determined. It is observed that with 
increase in the film thickness the mobility of the carriers increases 
whereas the carrier concentration decreases. The ase carrier 
mobility observed in the present study is 15cm? V's! on the 
er - (t=8000 A) with lowest carrier concentration of 7 
x10 : 


Effect of the thickness on the electrical properties of 

vacuum deposited AISb films is dealt with in this note. 

- Aluminium antimonide films were prepared by co- 
evaporation method as described by Francombe ef 

al.' All the deposition parameters such as substrate 
temperature, substrate, vacuum, etc. were held 

| _ constant to study the effect of the film thickness on 
“ioe their electrical properties. The film geometry used for 
_____ the Hall measurements was similar to that given by 
* Goswami and Ojha’. For the Hall measurements, a 
ae eee field of 9000 G as measured by a 


* ed _with digital multimeter by changing the 
ae | eee tes! titeetini. The film 


ps cn 


was applied. The Hall _ voltage was 


-alculatec pata eg a 


and the Hall mobility was calculated using the relation, 


h= R H.O 
where o is the measured conductivity. Fig. 3 shows the ; 
variation of y and u with thickness. It is observed that 
the mobility of the carriers increases with increase in 
the film thickness whereas the carrier concentration 
decreases. ; 


2000 4000 
t(A) 


Fig. 1—Hall constant (Ry) versus thickness of the 


a 
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NOTES 


The change in the electrical properties of vacuum- 
_ deposited thin film with varying thicknesses has been 
__ attributed mainly to the deterioration of the structural 
perfection. The electron microscopic observations 

show that the films prepared at the substrate 
temperature of 550°C were polycrystalline in nature. 
At lower thicknesses (t < 1000A) the crystallites 
would be smaller in size and are linked with the 
amorphous or polycrystalline material which offers 
high resistance. However, as film thickness increases, 
the size of the crystallites also increases and 
intercrystalline material offers less resistance and 
hence the conductivity of the film increases. 

It has been observed in the case of InAs single 
erystals* that the mobility dropped sharply when the 
carrier concentration varied from 1.7 x 10'° to2 x 
10'* cm ~>. Fig. 3 also shows, an increase in the carrier 
concentration with a decrease in the thickness which 
ultimately results in a decrease in the carrier mobility in 
the films. It may be mentioned that the carrier mobility 
in films drops with increasing carrier concentration 
more vigorously than in single crystals. This can be due 
to the presence of scattering effect in the films, whose 
contribution increases with decreasing film thickness. 
Similar observations have been reported by Wieder * in 
the case of InAs films, ie., a decrease in the film 
re “ aigapereniad us? ‘an increase in the dislocation 
- density which incre: s the dislocation scattering. It is 
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in thinner films. It has been reported® that the mobility 
decreases when the thickness of silicon samples is 
reduced by etching from 60 to 20 yu. 

Thus the increase in the carrier mobility with 
increase in film thickness can be attributed to a 
decrease in the carrier concentration, by a decrease in 
the intercrystalline barriers and dislocation scattering. 
The carrier concentration drops from 9 x 10!? em? 
at 2000 A thickness to 7 x 10!’ cm ~? at 8000 A. The 
mobility of the carriers increases with an increase in 
thickness and a maximum carrier mobility of 15 
cm?V ~'s ~! is observed during the present study on the 
thickest ‘fle (t = 8000 A) with the lowest cartier 
concentration of 7 x 10!’cm~3, 

The authors are indebted to Profs A R Patel and G 
K Shivakumar for valuable discussions. Thanks are 
also due to Drs S Chandra and C P Sharma of the 


National Physical Laboratory, New Delhi for useful 


discussions. One of the authors (AMB) is grateful to 
the University Grants Commission, New Delhi for 
the award of a Junior Resear) Fellowship. — 
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Thickness dependence of electrical resistivity of vacuum deposited 
selenium films (thickness, 500 to 4420 A) has been reported in 
vacuum, after an annealing cycle (30 to 150 C and back) and in 
residual vacuum after 24 hr. Unusual increase of electrical resistivity 
from 500 to 3400 A thickness has been explained in terms of localized 
surface states. 


The electrical conduction of thick amorphous 
selenium films has been reported to be space-charge 
limited'’? and in the lower temperature range (— 100 to 
20°C) it has been found to be polaronic’. Oxygen is 
known to increase the conductivity of pure selenium by 
many orders of magnitude in presence of even 50 ppm 


molecules of oxygen*’’. In the present note, we explain » 


the observed unusual increase in electrical resistivity of 
vacuum deposited selenium films with increasing 
thickness in vacuum, after an annealing cycle and then 
in residual vacuum after 24 hr. 

Specimens have been prepared by thermal 
evaporation of Se-powder of 99.5% purity onto 
cleaned glass substrates from a molybdenum boat in 3 
x 10° Torr vacuum at room temperature (30°C). The 
rates of deposition and film thickness have been 5 
A/sec and 500 to 4420 A respectively. Area of the films, 
0.2mm x 25mm has been obtained by shadowing the 
substrate by a copper wire during silver contact 
deposition. Leads from the silver contacts have been 


taken with the help of cleaned spring clips. Stable» 


values of electrical resistance have been measured just 
after deposition, after an annealing cycle (30 to 150°C 
and back) and in residual vacuum after 24 hr. Stable 
values of electrical resistance just after deposition have 
been attained in about 50 min. Rate of heating and 
_ cooling during annealing has been about 0.5°C/min. 


tWork done at Thin Films Laboratory, Physi Depa | 
University of Jodhpur, Jodhpur 342001 Kips secs i + 


$Present address: Department of Physics, K 
Saudi Arabia 
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(curve 2) and in residual vacuum after 24 hr (curve 3). is 


: of =2.2 1s Peas a 
Fig. | shows a typical behaviour of electrical x 10° oxygen molecules striking per unit area 


resistivity in the selenium films as a function of film. found that even 50 ppm of oxygen molecules reduce 


electrical resistivity of pure selenium (10"” ohm-cm) by 
-Tesistivities in the range 10’ to 10° ohm-cm, which 
ing Faisel University, i oxygen 


Log Resistivity (ohm-cm) 


2000 3000 4000 


Thickness (A) 


0 1000 


Fig. 1 — Electrical resistivity versus thickness [Curve: |, stable value 
after film deposition: 2, after an annealing cycle; and 3, after 24 hr in 
residual vacuum ] 


thickness. Electrical resistivity increases with thickness 
up to 3400 A and thereafter decreases. The initial 
increase in resistivity with film thickness is in 
contradiction to the thin film size-effect theory®’’. This 
unusual behaviour could be understood in terms of 
localized surface states, having energy distributed in 
the forbidden energy gap. These localized surface 
states exist due to chemisorption of oxygen molecules 
from environment, i.e. 2 x 10~° Torr partial pressure 
of oxyger® and non-crystalline structure of the 
specimens. These, in turn, cause statistical fluctuations 
in thinner region and form conducting channels, which 


are responsible for unusual behaviour below 3400 A 


thickness’. For films thicker than 3400 A, variation of 
resistivity is dominated by the size-effect on electron 
free path in the bulk and thus resistivity decreases 
thereafter. | 

No significant changes in the magnitudes of 
electrical resistivities of selenium films have been 
found for the following 3 conditions, viz. just after — 
vacuum deposition (curve 1), after an annealing cycle 


This could be due to adequate chemisorption of 


oxygen molecules during film deposition. In the 
‘Present study, the films have been deposited in an 


environment of 3 x 10~* Torr pressure (i.e. in presence 
of the substrate per second). LaCourse et al.* have 


six orders of magnitude. And our films have electrical _ 


confirms adequate chemisorption of oxygen. 
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A method for preparing multialkali photocathode starting with an 
antimony layer possessing island structure is described. The 
photocathode thus formed is found to be more sensitive in the near 
infrared region than that prepared by the conventional method. A 
reproducible photocathode with average yield of 160 #A/Im and 
having long wavelength threshold at 870 nm is obtained. 


Since the discovery of multialkali photocathode, 
Cs Na; K Sb, by Sommer’ in 1955, various methods of 
preparing the photocathode have been reported? es 
Efforts have been made in past to enhance the 
sensitivity of the photocathode in the near infrared 
region with a view to exploiting the predominance of 
infrared photons over the visible ones tn the night sky 
radiation’. Variation in the sequence of introduction 
of alkali vapours and different conditions of activation 
have been tried. Most of the methods have a common 
feature in that preparation of photocathode starts with 
the deposition of an Sb layer of controlled thickness 
which is then treated with different alkali vapours. 
Among the conventional photocathodes, Ag -O—Cs 
has the longest threshold extending up to 1.1 um 
(Ref. 5). This photocathode layer is known to consist of 
islands rather than being a smooth layer®. Takahashi” 
reported that the photocathodes prepared with initial 
antimony (Sb) layer in the form of a mosaic film rather 
than a smooth one had greater efficiency in the near 
infrared region of the spectrum. This note presents our 
experimental results on photocathode prepared from 
Sb layers having island structure. 

The photocathodes were deposited on a plain glass 
window attached to one end of a diode. The 
photocathodes were processed on an_ all-glass 
mercury-diffusion-pump-system which maintained a 
pressure of ~1x10~’ Torr during processing of 
photocathodes. An ion gauge, based on the 
phenomenon of thermal ionization of alkali atoms 
hitting an incandescent filament, was used for 
monitoring the partial pressure of alkalis during 
activation of the cathode. The details of construction 
of the diode and the experimental set-up have been 
reported elsewhere’. | 

After assembly, the diode was heated at 300°C 
overnight for 18h during which it was continuously 
pumped. When cooled subsequently to room 
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temperature the pressure in the system was registered 8 
x10~® Torr. Alkali generators were then degassed 
and the diode was again heated at 275°C for a further 
period of 18h and then cooled to room temperature. 

An Sb film having 30% light transmission was 
deposited on the glass substrate at room temperature. 

The diode was then heated to 250°C for about 15 min. 

The light transmission through the Sb layer increased, 
indicating thereby breaking-up of the smooth Sb layer 

into islands. Additional Sb was evaporated on the 
substrate till the light transmission of the film 
decreased by about 10%. The Sb film thus formed was 
treated with alkali vapours as follows: The diode 
temperature was raised to 180°C and potassium (K) 
was generated by resistance heating of the K- 
generator. White light from a calibrated tungsten lamp 

was utilized to monitor the photoyield of the 
photocathode: The variations of photocurrent and 
alkali ion current during the fabrication of 
photocathode are shown in Figs. | and 2 respectively. 
When the photosensitivity tended to. fall off the peak 
attained by K-activation, the diode temperature was 
raised to 200 C and sodium (Na) vapours were 
generated by heating the Na-generator. Sensitivity of 

the cathode then attained a new peak. As the sensitivity 
began to fall from the second peak, the photocathode 
substrate was air cooled to allow as much Na and K as 
possible to stick to the photocathode. Na- and K- 
generators were switched off when the photocathode 
current fell to 1/3 of the second peak. Sb was 
evaporated at that stage to increase the sensitivity. 
Subsequently, the photocathode was cooled by putting 

off the furnace. It was observed that sensitivity 
dropped to zero at room temperature. The tube was ; 
baked again and when the temperature was again 
raised to 100°C, potassium was generated. When the 
sensitivity attained a peak, Sb was evaporated on to 
the cathode till the sensitivity reduced to 1/3 of the 
peak. The sensitivity again attained a new peak. Again _ 
it was reduced to 1/3 of the peak by evaporating — 
additional Sb. K-Sb cycles were repeated in order to _ 
maximize the sensitivity. After obtaining the peak, the 
temperature of the tube was raised to 200°C and _ 
addition of Na was started. The evaporation of Na was 
continued till the photocathode sensitivity dropped to 
almost zero. The diode was then cooled to 160° 
K-Sb_ cycles were Tepeated to get 


Photoyield. K-generator was then switche: 
Cesium (Cs) vapours were then generated. Wh 
Sensitivity started falling from the peak, evaporation ¢ 
Sb at a controlled rate was started. \fter reaching 
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of the photocathode peaks at 570 nm, ies the 
wavelength threshold lies at 870 nm as compared v with — 
those reported by Sommer having peak at ee an ind ae 


“i * A a 


threshold at 820 nm (Ref. 11). ee a 
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A modified version of the famous Gell-Mann-Okubo mass 
formula for the pseudoscalar mesons has been derived from the (8,8) 
model of the broken chiral S U(3)@SU(3) symmetry. This modified 
formula agrees with the experimental data to within 0.4 %, which is 
much better compared to 6.3 °% for the original formula. 


The successes of the S U(3) symmetry in hadron physics 
are well known’ and need not be recounted here. It was 
realized later by Gell-Mann and others that the strong 
interactions are approximately invariant also under 
the larger group SU(3)@SU(3). For the symmetry 
breaking piece of the Hamiltonian, the most popular 
and successful model was the (3,3*)@(3*,3) repre- 
sentation?’> of SU(3)®@SU(3). 

However, it soon became apparent in the early 
seventies that the conventional (3,3*)@(3*,3) model 
cannot successfully account for certain experimental 
data, e.g. the meson-meson scattering lengths, meson- 
meson 0-terms, the amplitude for the y— 32 decay, etc. 
To remedy this situation, we propose the (8,8) 
representation of SU(3)®SU(3) as a_ promising 
alternative to the (3,3*)®(3*,3) model. The present 
authors have investigated the (8,8) model and carried 
out several calculations using it*” ’ with encouraging 
results. In this note, we exploit the (8,8) model to derive 
a modified version of the famous Gell-Mann-Okubo 
mass formula® for the pseudoscalar mesons. The 
formula derived here shows much better agreement 
with the experimental data than the original formula. 

Our approach is similar to that of Cicogna and 
coworkers” who analyzed in detail the consequences of 
the chiral symmetry breaking and discussed the results 


which follow simply from the (3,3*) ®(3*,3) model. For 


simplicity, we consider the case in which the field 
theory is described by a Lagrangian L(x) given by 


Ux) =Lolx) +49) (I) 


where Lo(x) is invariant under the given group G,d; are 
the real constants and ¢(x) are the basic local fields in 


: *De ERR of Matt tics, Regional Engineering College, 


terms of which the Lagrangian is constructed and 
which transform locally as a representation of the 
group G as 

P(X) gi; P(x) ae 
Here gi; is the matrix representation of the generator G* 


in the vector space spanned by the local fields ~ (x). 
Now we introduce the notation 


d=d,9; .. (3) 
where @; is the unit vector in the direction i, and A; are 
the vacuum expectation values of the local fields ¢,(x) 
and are given by 

A; = <0| @(x)|0> .. (4) 

With this notation, and using the functional method 
to treat the spontaneous breaking of symmetries in 
quantum field theory, it can be shown that the second- 
order Ward-identity can be written in the basis in which 
the mass matrix is diagonal in the form? 

m? EG", A\;=(G",d]; eee (5) 
where m? is the squared mass of the particle i. 

In the (8,8) model, the local fields g(x) are denoted 
by the symbols Z, and Z,. They are defined in the 
suggestive form of the products of SU(3)@SU(3) - 
currents“ as follows: , 

Zo =(1/2./2(V—A)§(V+ A) up 
Za=(/3/x/5)dap(V— A)g(V + A)y, 

Z,=(1//3) fapy(V—A)B(V + A)yy . 6) 
with «,8,y=1,2,....8; =0,1,2,3. 

In Eq. (6), the fields Zp and Z, belong, respectively, to 
the singlet and the even-parity octet in the S U(3) 
decomposition of 8@8, while Z, belongs to the odd- 
parity octet. The fields Z, correspond to the scalar 0* 
mesons, while the fields Z, correspond to the 
pesudoscalar 0° mesons. 

The strong interaction Hamiltonian H can be 
written as ; 
H=H,+H' lca: (7) 
where Hg is SU(3)@SU(3) invariant, while H’, which 
breaks this symmetry, can be written in terms of the 
fields given in Eq. (6), as 5 
H' = doz +dgZe ons (8) 

In Eq. (8), do and dg are the symmetry-breaking 
parameters. We have chosen the above form of 
symmetry breaking to satisfy the requirements of the 
conservation of isospin and hypercharge for the strong 
interactions. 

Now we utilize the Hamiltonian of Eq.(7) in the 
second-order Ward identity of Eq.(5) to derive the 


A=AiPj ’ 
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sSes ; scalar 
expressions for the squared masses of the same 
mesons*. By using the appropriate generators !or ee 
mesons 7, K and n and evaluating the commutators 0 


Eq. (5), we get 
/10do “4 2ds 
/10A¢ + 2Ag 
m2 2% \/10do —dg 

K </ 10A, == As 
m2 = ¥ 10d0— 24s .. 9) 

Ps Oka 2A 

Moreover, the expressions for the decay constants 
F, of the pseudoscalar mesons can be obtained easily 
by combining the hypothesis of the partial 
conservation of axial-vector currents (PCAC) with the 
Ward identity of Eq. (5). This gives: 
F, =C(./10A9 +2.) 
Fy =C(/ 1019 — As) 
F, =C(./10A9 — 2A) 
where 

C=(3,/3/10) 


Combining Eqs (9) and (10), we get the remarkable 
relation 


4F , Mx =3F ym; + Fam; 


= 


.. (10) 


oe 4a) 


which is the modified Gell-Mann-Okubo mass formula 
obtained from the (8,8) model. We note here that this 
relation is not new. It was obtained previously by 
several workers using various approaches. But it is for 
the first time that it has been derived from the (8,8) 
model and this shows once again that the (8,8) model 
may prove to be a promising alternative to the more 
popular (3,3*)® (3*,3) model. 

The original Gell-Mann-Okubo mass formula was 
derived in terms of linear masses from the S U(3) 
symmetry and proposed in the form _ 


4m, =3m,+m, _ (12) 


But the agreement obtained using this relation with the 
latest experimental data on the masses is very poor. It 
turns out that if we substitute the average experimental 
masses in Eq. (12), it is satisfied to within 10.6 %. 
The agreement can be improved to some extent if we 
use the squared masses in the place of linear masses, 
Then we get the relation 
4m? = 3m; +m; a (13) 
which is satisfied to within 6.3 vas However, the reason 
why the use of the squared masses improves the 
agreement is not understood jn depth. The only 
justification proposed for using m? instead of m 
(i=n,K.) was that in the relativistic field theories 
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of mesons (e.g. in the Klein-Gordon equation), the 


P 2 
meson mass always enters as Mm’. 
Now we confront our modified relation (11) with the 
experimental data and show that there is a marked 
improvement in the agreement. Using the accepted 


numerical values 


F,/F,=1.25 and F,,=0.95 m, ... (14) 
in the first two relations of Eq: (10), we get 

Jy =0.675 m, 

Ag= —0.152 m, xo LS 


With these numerical values, the first and the third 
relations of Eq. (10) give 


F, _ s/10d9—24s _ 435 


ale ii Rectan UT ce! 
Fr, ./10Ag+2Ag 


Substituting the numerical values of the meson 
decay constants from Egs (14) and (16) and the latest 
experimental data!° on their masses into Eq. (11), we 
find that our modified Gell-Mann-Okubo mass 
formula is satisfied to within 0.4%, which is really 
remarkable, since the chiral S U(3)@ S U(3) itself is only 
an approximate symmetry and its violation is expected 
to be of the order of 10-20%. 

Another satisfactory feature of our modified relation 
(11) is that in the limit of an exact SU(3) symmetry (i.e., 
F,~F, ~F,), it reduces to the original Gell-Mann- 
Okubo mass formula (13) with squared masses. But 
relation (13) suffers from the serious drawback that the 
squared meson masses were introduced in this formula 
on an ad hoc basis, merely to improve its agreement. In 
contrast to this, the squared masses appear in our 
modified relation (11) as a logical consequence of the 
(8,8) model and the second-order Ward identity, whic 
is more satisfying. 

We are thankful to the authorities and staff of the 
Meerut University Library for their help and 
cooperation. One of us (TNT) is grateful to the 
Ministry of Education, Government of India, for the 
financial assistance. tition ~” 
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The magnetic moments of charmed baryons are studied in the 
broken SU(8) framework, including the additional symmetry- 
breaking effects due to quark masses. The results so obtained differ 
from those of the conventional quark model incorporating quark- 
mass breaking in magnetic moment operator only. 


Magnetic moments have provided some very 
significant information in our search to understand the 
structure of the hadrons. The large moments of neutral 
baryons as well as anomalous magnetic moments of 
charged particles are clear evidence for important 
internal structure. The magnetic moments of baryons 
have been calculated in both the symmetry schemes! ~* 
and in the quark models® ~’ but now from the recent 
precise measurements’ it is clear that simple additive 
quark model considerations cannot explain these 
moments. Various attempts’ '° have been made to 
add symmetry-breaking contributions. De Rujula er 
al.**, ina mass broken quark model, assigning point- 
like Dirac magnetic moments to the quark which is 
taken to be proportional to the charge-mass ratio, and 
determining the quark mass ratio parameter from the 
baryon masses, successfully predictec (/)= 
—0.61 nm. But with the same quark mass ratio 
parameter the predictions for «(2*), u(=) and w(=°) 
are found to be much different. Identical results have 
also been obtained by Lipkin'* assuming that the 
quark mass difference is exactly equal to the baryon 
mass difference [(m,—m, = (Ma — M,)]. This implies 
that only the difference in quark masses does not 
provide the most general way to obtain symmetry- 
breaking contributions. Bag model'’ which is 
- 22 to be the manifestation of QCD and also 
| ; the relativistic effects does not improve the 
ion much. Then it was thought to add additional 
cata of the strong interactions in the 
tions. Fine structure corrections'*,such as 
tic effects, configuration mixing, isospin 
Z erin and pion and kaon cloud 
iO ? that one could expect in a constituent 

del, have been investigated, but 1% level of 
acy has not been achieved. Recently, Bohm ef 
LO toe calculated the octet baryon magnetic 
e ' yp seconc order perturbation theory, 
scoun pate effect of spin-dependent 


interaction arising from QCD, but failed to obtain an 
improved fit. 

Some workers?! have suggested different 
approaches of including the additional symmetry- 
breaking contributions. They have shown that the 
Situation improves if one brings in baryon mass 
dependence in the computation of baryon magnetic 
moments essentially by bringing in dimensional 
consideration in the symmetry predictions. The 
Situation appears to further improve?*?> when one 
considers the correction factor as given in terms of the 
Square root of the masses of the baryons. This situation 
can be understood because the physical mass reflects 
symmetry-breaking effect and therefore, such a scale 
factor correctly represents a symmetry-breaking 
property. In a recent paper, Tomozawa?° argued that 
such mass corrections may not be the right ones since 
symmetry-breaking effects presumably may be 
originating from the same source, (i.e. quark masses) 
and therefore, one should use the same parameters 
throughout. In calculating the octet baryon moments, 
he proposed?° a symmetry-breaking mechanism which 
keeps ground state SU(6) wavefunctions and 
introduces a mass scale factor which is a symmetric 
function of quark masses and has the dimension of 
inverse mass. 

The purpose of the present note is to study the 
magnetic moments of charmed baryons following the 
considerations of Tomozawa’° and also to make a 
comparison with the estimates obtained in other 
formalisms!?:75:7®, Since charmed baryons have been 
observed?’ and in the light of the possible mass 
spectrum some of the 4* charmed baryons are 
expected to be stable against strong interactions, such 
studies may help in deciding the pattern of symmetry- 
breaking mechanism. 

The magnetic moment operator'* ~'® pis expressed 
as 


p= 0 .. (D 
where o=(3-4-3+) pe 6. 
with A=, x= Ft and bone rey | 


and m,, Mg, m, and m, are the masses of respective 
quarks. The magnetic moment of charmed baryons 
can be calculated by taking the expectation value of the 
operator (Eq. (1)] between spin unitary spin 
wavefunctions’ and are given in the second column of 
the Table 1. For the additional symmetry-breaking 
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Table 1—Charmed Baryon Magnetic Moments 
Particle SU(8) matrix Mass factor 
element of 


Eq. (1) 
(4-9) 
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term coming from strong interaonene we: 
mass scale factors 


l 
(A) ge 


(B) +h tas 
where ™,, 2, M3 are the masses of the quar 
constitute the baryons. The mass scale factors 
and (5)] entering in the expression for the n 
moments of various baryons can easily be ca 
and are given in third and fourth columns of ’ 
The parameters 4 and x can be determined 
4p) = 2.79 nm, 4(n) = —1.91 nmand Ny: = 
nm as inputs and these come out to be sia 
=().726 (Ref. 26). For the parameter y, as 
charmed baryon whose magnetic moment 
we following Rujula eral.'* etgsee 
the known baryong masses in the € 


2My++ My,— 3M x. 
~ 0.23 3 


y= 


Once the ooo ie se 


NOTES 


results are comparable. However, it will presumably 

take some time before the magnetic moments of 

charmed hadrons can be measured. Nevertheless, we 

present our results now because, there may be other References 

calculations more accessible to experimental tests in —_1 Coleman § & Glashow S L, Phys Rev Let (USA), 6(1961) 423. 

which magnetic moments are needed. Also, the decay 2 BegMA, Lee B W & Pais A, Phys Rev Leti(USA), 13(1964) 514. 
ca Tates of some of the charmed baryons willdependupon 3 Choudhary A L & Joshi V, Phys Rev D (USA), 13 (1976) 3115. 
wee i 
5 


-The authors are thankful to Prof. B S Bhargava for 
providing the necessary facilities in the Department. 


: transition moments which may be of immediate Verma R C & Khanna M P, Pramana (India), 8 (1977) 462. 


: ' Rubinstein H R, Scheck F & Socolow R H, Phys Rev(USA), 154 
In the following, the relations among different ratios (1967) 1608. 


of the magnetic moments of the baryons of the same Franklin J, Phys Rev (USA), 172 (1968) 1867; 182 (1969) 1607. 
a ‘multiplet d Singh L P, Phys Rev D (USA), 16 (1977) 158. 
re are discussed. However, if we neglect the 


E Particle data group, Rex Mod Phys (USA), 52 (1980) 1. 
__ isospin breaking, i.e. take A=1, the following three 9 Kamal AN, Phys Rev D (USA), 18 (1978) 3512. 


s are obtained: Singh C P, Verma RC & Khanna M P, Pramana (India), 13(1979) 
261. 
Pandit P N, Singh C P & Khanna M P, Acta Phys Austr( Austria), 
§3 (1981) 211. 
Verma R C, Phys Rev D (USA), 22 (1980) 1156. 


Rujula A De, Georgi H & Glashow § L, Phys Rev D(USA), 12 
(1975) 147. 
Lipkin H J, Phys Rev Lett (USA), 41 (1978) 1629. 
Franklin J, Phys Rev D (USA), 20 (1979) 1742: Phys Rev Lett 
(USA), 45 (1980) 1607. 
= ee ae : Lichtenberg D B, Phys Rev D (USA), 15 (1977) 345; Z Phys s 
os < az ) (Germany), 7 (1981) 143. ‘ 


x ie): = (=) es 17 Allen E, Phys Lett B (Netherlands), $7 (1975) 263: Hackman R H 
ee | et al., Phys Rev D (USA), 18 (1978) 2537. : 
18 Isgur N & Karl G, Phys Rev D (USA), 21 (1980) 3175. — 
19 Theberge S & Thomas A W, Phys Rev D (USA), 25 (1 
Gozizaliz P et al., Kaon cloud and ss 
~moments, D Ph Til 19, 1982. 
20. Bohm M, Hwerta R & Zepeda A. 


_ Universi preprint 1980. es 
Meneed Rev eerie I m 


Indian Journal of Pure & Applied Physics 
Vol 21, July 1983, pp. 430-431 


EPR of Cu?* in MgSeO,.6H20 
Single Crystals 


V K JAIN 
Department of Physics, M D University, Rohtak 124001 
Received 27 November 1982; revised received 27 January 1983 


The EPR spectrum of Cu?* (S=1/2, /=3/2) in MgSeO,.6H,0 
consists of two overlapping angle-dependent four-line hyperfine 
patterns arising from two differently but magnetically equivalent 
Cu?* complexes located at a specific site. The angular variation of 
the spectrum indicates that Cu2* substitutes for Mg?” at Y-ion sites 
only. 


In the study of complex crystal structures in which the 
isomorphic substitution of an impurity ion can occur 
at a number of different lattice sites, the general 
methods of crystal chemistry are often inadequate 
when it comes to predicting which crystallographic 
position must be occupied by the ion. Many question 
pertaining to the mtroduction and the state of impurity 
ions can be elucidated from the electron paramagnetic 
resonance (EPR) spectra. The X-ray analysis of 
MgSeO,.6H,O single crystals has shown that the 
crystal structure consists of two crystallographically 
different Mg** ion sites’. Subramanian et al.? and Jain 
and Venkateswarlu® have studied respectively the EPR 
of Co** and Mn’* in MgSeO4.6H,0O single crystals 
and shown that impurity ions occupy both the sites. 
The present study deals with the results obtained in the 
EPR of Cu** incorporated into MgSeO,.6H,O single 
crystals. 

Magnesium selenate hexahydrate (M gSeO,.6H,O) 
crystallizes in monoclinic symmetry with space group 
C(2/c) — C$,!. The lattice constants are a=10.36 A, b 
=7.38A, c=25.1 A and B=98.1°. The unit cell 
contains eight formula units. Out of eight Mg?* 
Present in the unit cell, four have the same 
crystallographic surrounding having point symmetry 2 
(X-ions). The other four have point symmetry I (Y- 
ions). The Y-ions form two Pairs which are equivalent 
through reflection in the ac plane. The Mg?* at both 
the sites are surrounded by a distorted octahedron of 
water molecules. 

Crystals of MgSeO,.6H,O doped with Cy2+ were 
grown at room temperature (~ 298 K) by slow 
evaporation of saturated aqueous solution. The Cy2* 
was introduced into the host lattice by adding a small 
amount (1% by weight) of cupric sulphate. EPR 
experiments were performed on a JEOL FE-3X 
homodyne X-band spectrometer equipped with 
TEo;;-cylinderical cavity and 100 kHz field 
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by the space group symmetry operations. For X-ions, 


modulation. As a reference for magnetic field strength, 
the resonance line of DPPH with g = 2.0036 was used. 
The crystals were mounted on quartz rods. The 
angular variation studies were done using JES-UCR- 
2X sample angular rotating device. 

* The EPR spectrum of Cu?” (S=1/2, 1 =3/2) in 
MgSeO,.6H,O consists of two overlapping angle- 
dependent four-line hyperfine patterns (Fig. 1) arising 
from two differently but magnetically equivalent Cu*~ 
complexes located at a specific site. When the magnetic 
field is along the b axis or in the ac’ plane (c’ is at 90 toa 
axis in the ac plane) only a single set of four lines was 
observed indicating that none of the principal axes lie 
in ac’ plane or along the b axis. The angular variation of 
the spectra in ub, bc’ and ca plane shows orthorhombic 
symmetry. The linewidth of the spectrum along the = 
axis is ~2.0 mT. The spectra at 298 K have been 
analysed using the spin-Hamiltonian 


Sf= tp S.2.B4+8S.A.1 


where the symbols have their usual meaning. The 
values of the spin-Hamiltonian parameters are: A- 
=105.84+1.0, A,=474+1.0, A, =1664+2.0, g. 
= 2.400+ 0.001. g,=2.109+0.002, g, = 2.080 
+ 0.002, 0=80 +1°, p= —20°+2° where the A 
values are in units of 10° “ecm '" 0 is the angle between 
the = axis and the b axis and ¢@ is the angle between the a 
axis and the projection of z on the ac’ plane. The values 
of the components of the g-factor indicate that Cu2* in 
the lattice is surrounded by six water molecules‘. 

If Cu** enters substitutionally in MgSeO,.6H,O, 
the sites available to it are those of Mg?* and Se®*. As 
the valence states of Cu?* and Mg** do not differ, 
Cu’* would prefer to occupy Mg?” sites. The crystal 
structure of MgSeO,.6H,O requires that Cu?* 
oriented arbitrarily in both X- and Y-ions site have 
four symmetry-related equivalent orientations given 


* 


all are magnetically equivalent for all orientations of 
the crystal. All Y-ions are chemically equivalent but __ 
magnetically equivalent in pairs, showing two sets of 
quartets for any arbitrary orientation. The angular 
variation of the spectrum indicates that Cu2* 


° ee I OOmT, . — 
Fig. |—EPR spectrum of Cu2* in Gisele 
sh ace ame, ome 
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substitutes for Mg** at Y-ion sites only. Because if it 
only substitutes for X-ion sites, the spectrum should 
consist of four lines for all orientations of the crystal as 
required by the space group symmetry operations. 
Nevertheless, if Cu?” substitutes for both X- and Y- 
ions sites, one would have observed two sets of four 
lines for magnetic field in the ac’ plane or along the b 
axis. Further, a comparison of the value of 6 and 
obtained experimentally with those calculated from 
the data of Zalkin et al. for different Mg—H,O bonds 
indicates that z axis is nearly along the longest Mg (Y- 
ion)\—H,O bond, ie. Mg—H,O(7) (0=79.0, p= 
— 14.36"). Thus Cu?* substitutes for Y-sites of Mg”* 
to form [Cu(H,O),]?* complexes with z axis nearly 
along the longest Mg—H,O bond direction. On the 
other hand, EPR of Co** and Mn’* in 
MgSeO, .6H,O indicates that impurity ions enter the 
lattice substitutionally at both X- and Y-ions sites". 

The [Cu(H,O),]?* ion in the MgSeO,.6H,O 
belongs to the point group 1(C,), so that the metal part 
of the ground state wavefunction in the absence of 
spin-orbit coupling consists of an admixture ofall five d 
orbitals. In fact, the ligand field is composed of a 
dominant D,, component with an additional 
perturbation due to the inequivalent bond length along 
the x and y axes which lowers the point group to Dp 
and a further lowering of the symmetry to C; because of 
the non-planarity of the atoms in each Cu—H,0 unit. 
In the D,, point group, the metal part of the ground 
state is de_,2 and lowering of the symmetry to D5) 
causes this to become contaminated with d3.--. The 
further lowering of the symmetry to C; is accompanied 
by a change in the nature of the interaction of the water 
molecules with the metal d,,, dx-, and d,- orbitals and 
consequent participation of these in the ground state. 
However, participation of d,,, d,., d,. orbitals have 
little effect on the g values. Thus neglecting the 
participation of d,,, d,. and d,. in the ground state 
wavefunction, the metal part of which can be written° 


y = a|d.2-,2)> = b\d32— > 
where a2 + b2 =1. The equations for the g shifts are 
those for a complex of D,, symmetry”’° 


Ag, = — 2(A/E,.)k{(a — ./3b)? 

Ag, = —2(A/E,.)k3(a + ./3b)’ 

Ag, = — 8(A/E,,,)k3a? 

where EF; represents the energy of the excited state in 


which the unpaired electron cccupies the orbital 
denoted by the subscript; A is the spin-orbit coupling 
constant (—828 cm ' for Cu’*); and k,, k,, k3 are 
orbital reduction parameters’. The net effect of the 
admixture of the d3--_,: into the ground state is to 
produce a new wavefunction of the form W |cx? 
+e)? + fz?) with unequal lobes along x, ), and z 


where c=b + \/3.a,e=b— /3.a, f= —2b.° 

As the electronic spectrum of Cu** in 
MgSeO,.6H,O is unknown the value of E,., E,. and 
E,, have been taken as those for [Cu(H,O),]** 
complexes in other hydrated crystals®, i.e. E,. = E,. = 
~ 12,500 cm™',: and -E,,=~10,000cm™'. The 
parameters c, e and f were calculated with the 
restriction k, =k, =k, and k; =k). The parameters c, 
e and f calculated from the molecular g values are c 
= 1.776,e = — 1.684, and f = — 0.092, k,? = 0.602, k ,” 
= 0,691: 
This work was supported by CSIR, New Delhi under 
project 3(502)/81/EMR-II. The author wishes to thank 
Prof. T S Jaseja for his interest and encouragement. 
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Twyman Green interferometer with curved mirrors is described. It 
is found that the end mirrors in the Twyman Green interferometer 
need not be necessarily plane but must be identical in figure. 


Twyman Green interferometer’ is used for testing the 
optical homogeneity of prisms, lenses and parallel 
plates in optical testing laboratories. It is essentially a 
modification of the Michelson interferometer where 
the extended source of light is replaced by a point 
source at the focus of a well corrected lens and 
consequently a collimated beam of light is used in the 
interferometer. The fringes obtained are therefore 
straight and equidistant for the case when the real 
mirror and the virtual mirror are making a small wedge 
angle. Fringe-free field is obtained for the case when 
the real and virtual mirrors are parallel to each other. 
When one of the mirrors in it is replaced by a slightly 
curved (convex or concave) spherical mirror, a set of 
circular fringes is obtained in the field of view due to 
the interference of the plane and curved wavefronts. 
The number of fringes depends on the amount of 
curvature present in the mirror. This type of 
arrangement is not suitable for testing the parallel 
plates and prisms. It is interesting to see the nature of 
the fringe pattern when both the plane mirrors are 
replaced by slightly curved (convex or concave) 
spherical mirrors which are identical in shape. Such 
type of arrangement is shown schematically in Fig.1. 
The mirrors are chosen to be about 2.5 wavelengths 
convex and placed at equal distances from the beam 
dividing surface. The reflected wavefronts from both 
the mirrors are identical in nature. When the 
interfering beams are parallel to each other, a uniform 
illumination Is obtained Over the aperture. In order to 
obtain straight fringes a slight tilt between the beams 
can be introduced. Fig.2 shows the photograph of 
nearly straight fringes obtained with the arrangement 
shown in Fig. 1. The result shows that it is not necessary 
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to make the mirrors perfectly plane but efforts must be 
made to make both the figures identical. 

The maximum amount of curvature that can be 
tolerated on the mirrors can be derived from Figs.(3a) 
and (3b). Fig.3a shows the situation when the two 
mirrors are not located at equal distances from the 
beam divider. This happens while testing the ruby 
rods, thick parallel plates, etc. The image of the mirror 
M, in the beam divider is separated by an amount / 
from the mirror My. The reflected wavefronts have 
their foci at F, and F,. The interference takes place 
between two wavefronts, whose foci are displaced 
longitudinally and therefore the illumination will not 
be uniform over the whole aperture. The variation of 
optical path difference (0) over the aperture is given by 
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